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This work contains a systematic treatment of single particle Synchrotron Radiation and some ap- 
plication to realistic beams with given cross section area, divergence and energy spread. Standard 
theory relies on several approximations whose applicability limits and accuracy are often forgotten. 
We begin remarking that on the one hand, a paraxial approximation can always be applied without 
loss of generality and with ultra relativistic accuracy. On the other hand, dominance of the accel- 
eration field over the velocity part in the Lienard-Wiechert expressions is not always granted and 
constitutes a separate assumption, whose applicability is discussed. Treating Synchrotron Radia- 
tion in paraxial approximation we derive the equation for the slow varying envelope function of the 
Fourier components of the electric field vector. Calculations of Synchrotron Radiation properties 
performed by others showed that the phase of the Fourier components of the electric field vector 
differs from the phase of a virtual point source. In this paper we present a systematic, analytical de- 
scription of this phase shift, calculating amplitude and phase of electric field from bending magnets, 
short magnets, two bending magnet system separated by a straight section (edge radiation) and 
undulator devices. We pay particular attention to region of applicability and accuracy of approx- 
imations used. Finally, taking advantage of results of analytical calculation presented in reduced 
form we analyze various features of radiation from a complex insertion device (set of two undulators 
with a focusing triplet in between) accounting for the influence of energy spread and electron beam 
emittance. 

PACS numbers: 41.60.Ap, 41.60.-m, 41.20.-q 



I. INTRODUCTION 

About sixty years have passed since the first, pioneer- 
ing works on Synchrotron Radiation have been published 
(see |l| and references therein). The way scientists con- 
sider this phenomenon has drastically changed during 
this period. At first, Synchrotron Radiation was regarded 
only as a detrimental factor, a limitation on the maximal 
particle energy attainable with accelerators (see, for in- 
stance, |2|). Nowadays it is an outstanding research tool 
allowing continuous advancement of both fundamental 
and applied sciences and it is used worldwide by tens 
of thousands of scientists from many different disciplines 
like physics, chemistry, material sciences and structural 
biology. 

The properties of Synchrotron Radiation can be de- 
rived by applying the methods of classical electrodynam- 
ics to the motion of relativistic electrons (or positrons) in 
magnetic structures like bending magnets or undulators. 
The theory of Synchrotron Radiation based on the works 
by Iwanenko, Pomeranchuk and Schwinger has been de- 
scribed in detail both in convenient summaries as well 
as in textbooks 0, S Ifl S S IE El • 

On the one hand, these texts contain derivation of fun- 
damental equations describing spectral-angular proper- 
ties of Synchrotron Radiation in the far-zone for single 
electrons. A usual starting point for quantitative descrip- 
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tions of these properties is provided by the formulas for 
the Lienard-Wiechert fields Q|. 

On the other hand, in relation with the construction 
of beam lines at Synchrotron Radiation facilities it is im- 
portant to understand the modifications to single particle 
treatment which is due to realistic beams with given cross 
section area, divergence and energy spread. 

This raises several questions. First, discussing finite 
electron beam emittance effects, an obvious remark is 
that, in the case of second generation light sources, the 
emittance of the electron beam was much larger than the 
radiation wavelength of interest. As a result, both pho- 
tons and electrons phase space had the same hamiltonian 
structure, and properties of radiation could be treated on 
the basis of geometrical optics. In recent years though, 
more and more specialized and optimized Synchrotron 
Radiation sources have started operation. During de- 
sign and optimization of these sources much emphasis 
has been directed at reducing the transverse phase space 
area of the electrons. Among the most exciting proper- 
ties of today's third generation sources is the tiny vertical 
emittance which is an order of magnitude smaller than 
1 A (for design parameters of an up-to-date source see, 
for instance, [HI). In this case geometrical optics cannot 
be applied any more and even basic properties like angu- 
lar distribution of intensity at fixed frequency, or spectral 
intensity distribution at fixed observation angle should be 
calculated applying full electromagnetic equations. 

A second question arises considering the features of un- 
dulator insertion devices. In older sources, typical length 
of insertion devices was about 1 2 m, while the dis- 
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tance to the user station was of order 20 -j- 30 m. In this 
situation the asymptotic formulas for the far zone could 
be applied and the question of their applicability region 
was of theoretical interest only. Nowadays installation of 
20 meter-long undulators has been planned at PETRA 
III (see 0) and distances between source and observer 
of several tens of meters have been proposed. A similar 
device (a 25 m long segmented undulator with no focus- 
ing elements between the segments) has been installed 
and operates at SPring-8 (see 0). In this situation ap- 
plicability region of far-field formulas and correction for 
near field effects are of actual practical importance, not 
to mention the interest for far infrared edge radiation. 

Another problem is related with the increasing com- 
plexity of insertion devices. For instance at PETRA III 
(see again 11]), installation of two undulators segments 
separated by a focusing triplet in between has been pro- 
posed to obtain small vertical betatron functions. Com- 
putation of radiation characteristics from these novel de- 
vices constitute a rather challenging problem. 

Furthermore, a question which is somewhat related to 
all the previous ones is linked with the dramatic increase 
of brilliance achieved in third generation light sources 
with respect to older designs, which has triggered a 
number of new techniques and experiments unthinkable 
before. Among the most exciting properties of today 
third generation facilities is the high flux of coherent x- 
rays provided. The availability of intense coherent x-ray 
beams has fostered the development of new coherence 
based techniques like fluctuation correlation d yna mics 
[3 13 13 13 13 Dl, phase imaging [3, M, M El, 
and coherent x-ray diffraction (CXD) |23l I24. l25l l'26| . 
In all these fields, the understanding of the evolution of 
transverse coherence properties along the beam line is 
of uttermost importance: in particular, both the beam 
size and divergence should be taken into account in the 
evolution model. 

As one can see these questions address quite different 
physical phenomena; yet they all belong to the field of 
Synchrotron Radiation. This number of different issues 
stems from the fact that practical applications make use 
of a very wide range of radiation wavelengths, from 0.1 A 
to 100 /im, which span over seven orders of magnitude. It 
is no surprise that very different problems arise when the 
wavelength of interest is tuned to such different regions of 
the electromagnetic spectrum. Considering, for instance, 
bending magnet radiation, it is clear that near field effects 
will be of practical importance for FIR (Far Infra-Red) 
applications but hardly for X-ray radiation, because of its 
much shorter radiation formation length. On the other 
hand, the influence of the omittance on X-ray radiation 
will be important because the magnitude of the omittance 
is comparable with the radiation wavelength, but it will 
be negligible when it comes to the characterization of FIR 
radiation properties, since in this case the wavelength of 
interest is way larger than the emittance. 

As a result, all the matters mentioned above are of 
great practical importance in some particular region of 



the spectrum, although they are systematically neglected 
in Synchrotron Radiation textbooks and reviews. Reac- 
tion of synchrotron-radiation-users communities to this 
set of problems was the development of computer codes 
capable to calculate radiation properties from realistic se- 
tups starting from first principles [23, HI . At first glance 
this settles all issues since, at least from a practical view- 
point, users have the possibility to calculate all the prop- 
erties they need. 

Yet, computer codes can calculate properties for a 
given set of parameters, but can hardly improve physi- 
cal understanding, which is particularly important in the 
stage of planning experiments. Understanding of correct 
approximations and their region of applicability can sim- 
plify many tasks a lot, including practical and non-trivial 
ones. The use of dimensional analysis in an analytical 
framework is of uttermost importance to this goal. By 
that, dimensionless quantities of physical interest can be 
easily identified. In particular, dimensionless parameters 
much smaller (or much larger) than unity always corre- 
spond to simplifications of fundamental equations: most 
physical theories originate from recognizing and taking 
advantage of such small (or large) quantities. 

In contrast, relying only on computational power to 
solve physical problems, besides being inelegant, would 
automatically neglect the presence of small (or large) pa- 
rameters and simplifications related to them. The lack, 
or poor understanding of approximations involved in a 
certain theory will result, in its turn, in poor understand- 
ing of the mechanisms involved in physical phenomena 
which would undoubtedly result in gross mistakes. More- 
over, just a few codes are available, which can be mod- 
ified to account for particular situations (like coherence 
properties or radiation characteristics from a complex in- 
sertion device for instance) by a few experts only, while 
there are about twenty facilities each with its own set of 
particular setups and situation to be accounted for. 

We propose to overcome all these difficulties by devel- 
oping a theory of Synchrotron Radiation which accounts 
for all complexities analytically. Analytical approach will 
help to understand physics of real Synchrotron Radiation 
sources, thus filling the gap between idealized textbook 
analysis and actual situations. We will make explicit use 
of small (or large) parameters involved in the description 
of cases of practical interest. In this way we will reduce 
problems as much as possible, to a level where physical 
insight is easy, so that simple computer codes, capable to 
describe the situation under study, can be developed also 
by non-expert programmers. In this process we will put 
particular care in the specification of the region of appli- 
cability and of the accuracy of the approximations made, 
which is usually neglected in analytical calculations. 

In the next Section we discuss methods for computa- 
tion of Synchrotron Radiation from a single particle on 
a generic trajectory. We start reviewing the two algo- 
rithms used today, the first based on Lienard-Wiechert 
fields, the second on Lienard-Wiechert potentials. We 
propose to start calculations from the very beginning 
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by solving paraxial Maxwell equations for a given har- 
monic of the field. This is, from a logical and educational 
viewpoint, the most natural way to perform calculations; 
from a computational viewpoint it is also the simplest 
way. The use of a paraxial approximation is justified 
by the features of an ultra relativistic system character- 
ized by a Lorentz factor 3> 1: radiation formation 
length is, then, much longer than the wavelength, and 
the radiation is distributed within a cone with opening 
angle much smaller than unity. Maxwell equation for a 
given harmonic of the field E(r,uj), which has elliptic 
characteristic, can be then transformed into a parabolic 
equation. We solve this equation with the help of an ap- 
propriate Green's function thus coming, in Section [HJ to 
a very generic expression for E. We end the Section by 
addressing the obvious question of the relation between 
our method, and the two algorithm introduced at the 
beginning. 

Subsequently we apply our expression to recover 
some well-known and less well-known properties of Syn- 
chrotron Radiation from bending magnets and undula- 
tors, in Section IllTl and Section llVl respectively. This is 
far from being a mere repetition of already known results. 
In fact our method will prove superior when it comes to 
physical understanding and determination of applicabil- 
ity region and accuracy of approximations: this has, of 
course, important application in estimation of practical 
quantities, like field intensities, with their accuracy. We 
will also discuss in detail the methodological issues about 
definition of fields and intensities from single devices. As 
an example of application of our new understanding we 
will show, in Section how radiation characteristics 
from a complex setup can be discussed in very simple 
terms. Finally, in Section IVII we come to conclusions. 



II. 



METHOD 



A. Review of known methods 

There are two basic methods which are used to calcu- 
late Synchrotron Radiation characteristics. 

The first [1 II S, 0, S, II 03 iS b&Sed ° n Lienard " 
Wiechert fields. It is the better known of the two and it 
is used in a very widespread way in textbooks. Lienard- 
Wiechert fields can be customary separated in a velocity 
and an acceleration term. Usually, the acceleration part 
alone is analyzed in Fourier components of the electric 
field vector: 



E(r„,w) = — 



dt'- 



(n - /3) x f3 



cJ-oo \r -r'(t')\(l-n-PY 
x.p/tofr+ l-- (1) 



the particle position, r D is the observer position Pit') 
and P(t') are, respectively, the particle velocity, v(t'), 
normalized to the speed of light c and its derivative, all 
calculated at time t' . It is clear that use of Eq. implies 
that the so-called velocity part of the field is neglected in 
the computation. 

The second method US, HI El HO] is based on Lienard- 
Wiechert potentials. As an alternative to the use of 
Lienard- Wiechert fields as a starting point for computa- 
tions based on first principles, a few authors of scientific 
research papers start with the Lienard- Wiechert poten- 
tials (A(r ,t), cj){r ,t)). These can be decomposed in 
their harmonics (A(r , w), <p(r a , u>)) which can be subse- 
quently used to calculate E{r ,uu) according to 



E(r„,u) = 
P — n ic 



iuje 
c 



dt' 



n 



|r„ - r'(t')\ a |r„ - r'(t')|= 
xexpW t - + l"°-'"'"'|l. (2) 



where (— e) is the electron charge, n(i') is the unit vec- 
tor from the particle to the observer position, r'(t') is 



We failed to find textbooks using this second method for 
educational purposes: although some authors (see [13 
for instance) start their derivation with harmonic analy- 
sis of potentials, when it comes to actual calculation of 
Synchrotron Radiation they go back to the first method 
based on Lienard- Wiechert fields. 

Questions arise about the relation between Eq. 
and Eq. ©. First, in contrast to Eq. JTJ, Eq. © 
is exact, in the sense that the velocity field term is not 
neglected. It is then natural to ask what are the condi- 
tions for which the velocity field term in Eq. Q can be 
dropped, and with what accuracy this can be done. In 
fact, when approximated expressions like Eq. Q and its 
specialization to particular magnetic systems (like un- 
dulators and bending magnets) are found, special care 
should be taken in specifying their region of applicability 
and accuracy; yet we failed to find, through references 
P. 13. Ifj HI IsL la. HfT| quantitative specification of accuracy 
of results obtained under given approximations. Second, 
Eq. and Eq. J5J are used both in the far and in 
the near field region and simplified under the paraxial 
approximation: it is important to understand relations 
between the paraxial approximation and the far or near 
field assumptions as well as region of applicability and 
accuracy of results. Third, Eq. is often used to- 
gether with integration by parts and neglecting the edge 
terms after integration: in this case it is interesting to 
ask what are the conditions which allow the edge terms 
to be dropped. Moreover, as a remark to both methods, 
we should note that, in general, one needs to know the 
entire history of the electron from t' = — oo to t' = oo 
since integration in Eq. and Eq. J5J is performed be- 
tween these limits. This statement should be interpreted, 
physically, depending on the situation under study: inte- 
gration should in fact be performed from and up to times 
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when the electron does not contribute to the field any- 
more. This does not present any problem in the case, for 
instance, of a circular motion because the particle tra- 
jectory is limited in space. In practical situation though, 
one deals with radiation from a beamline which is a com- 
position of insertion devices, like undulators and bending 
magnets, and straight sections. The requirement of in- 
tegration over the entire history of the electron poses, 
then, the following methodological problem: one has to 
give a meaning to the concept of radiation from a given 
insertion device, independently of the trajectory followed 
by the particle before and after the device. This problem 
will be discussed in more detail as we go through our pa- 
per. In general one can break up the electron trajectory 
in several parts, corresponding to physical devices like 
undulators, drifts and bending magnets, thus calculating 
the integral in Eq. or Eq. J5J along a finite time 
interval which corresponds to a particular device. In do- 
ing this, it is assumed that the observer is located far 
away enough from the device, so that the velocity part 
of the field is negligible with respect to the acceleration 
part. In this way, again, dropping velocity field part in 
Eq. m is justified and the results from Eq. and Eq. 
(j2J coincide. It is understood that different contributions 
from different devices must be accounted for with the cor- 
rect relative phase relation to get the total E(r,u>) as a 
sum of the separate contributions from single devices. In 
this way concepts like radiation field from a undulator 
or radiation field from a bending magnet have a well de- 
fined, practically useful meaning as partial contribution 
to the total field at the observation point. The concept 
of intensity is, instead, subtler, since it involves calcula- 
tion of | E(r,cj) \ 2 . In this case then, also interference 
terms between different devices should be accounted for, 
in the most generic case. Only when these are not im- 
portant the total spectral intensity from a given beamline 
can be broken up in the sum of separate terms from the 
single devices, and radiation intensity from a undulator 
or radiation intensity from a bending magnet are, then, 
well-defined. 



Finally, Eq. is based on the integration of Maxwell 
equation in time and in the subsequent harmonic anal- 
ysis. This is quite involved from a logical viewpoint: it 
is much simpler to start with the equations for a given 
harmonic u) so to obtain directly E after integration of 
Maxwell equations. A similar comment can be made re- 
garding the way Eq. J2J was obtained: first the potentials 
were found, then their Fourier transform was taken. In 
both cases one is bound to solve the full Maxwell equation 
even though the ultra relativistic nature of the systems 
considered in practice allows systematic application of 
a much more convenient paraxial approximation. Note 
that, usually, paraxial approximation is applied, but only 
at a later stage, which often makes analytical derivations 
more involved, and region of applicability and accuracy 
of results difficult to specify. 



B. Derivation using paraxial Green's function 

A system of electromagnetic sources can be conve- 
niently described by its charge density p(r, t) and current 
density j(r, t). In this paper we will be concerned about 
a single electron so that, using the Dirac delta distribu- 
tion, we can write 



and 



p(r,t) = -eS(r - r'(t)) 



j{r,t) = -ev{t)S{r-r'(t)) , 



(3) 



(4) 



where r'(t) and v(t) are, respectively, the position and 
the velocity of the particle at a given time t in a fixed ref- 
erence frame. We will assume that the particle trajectory 
is given a priori, and that Maxwell equations in vacuum 
allow a good description of the source field. In this case 
the magnetic field B can be easily recovered when the 
electric field E is known simply remembering 



B(r ,t) = -nxE 

c 



(•5) 



where n is the unit vector pointing from the retarded 
position of the particle, r'{t') to the position T"o at time 
t = t' + r'(t')/c. Therefore in this paper we will concen- 
trate on the electric field E(ro,t). 

Trivial manipulations of Maxwell equations in cgs sys- 
tem: 



and 



„ „ ldB 

VxE = — — 

c at 



„ 4vr . 1 dE 

c c at 



give the well-known equation for the electric field: 



c 2 V x(Vx£) = - 
With the help of the identity 



d 2 E dj 

47T — 

dt 2 dt 



(6) 



(7) 



(8) 



Vx(Vx£) = V(V-£)-V 2 £ (9) 

and Poisson equation 

V • E = Aup (10) 
we obtain the inhomogencous wave equation for E 



o , d 2 E dj 

c 2 V 2 E - = 4ttc 2 Vp + 4tt-^ 
ot at 



(11) 
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Eq. (|llfl is a partial differential equation of hyperbolic 
type. We introduce the Fourier transform f(to) of a quan- 
tity f(t) as follows: 



so that 



(12) 



f(u) = dtf(t)e wt , 



f(t) = — I dw/H 



(13) 



Using the representation in Eq (|13|) for the quantities in 
Eq. ((TTJ we get 



c 2 X7 2 E 



-u?E 



4irc Vp — iwiojj 



(14) 



Eq. (|14|l is the well-known Helmholtz equation which 
has elliptic characteristic. Upon introduction of a fixed 
cartesian reference system (x,y,z), it is always possible 
to define a quantity E such that 



Ee luz/c . 



(15) 



However this definition is useful only in the case when 
E varies slowly along z with respect to the length A = 
2irc/u}: only in that case, in fact, Eq. (|15f) is a useful 
factorization of E as the product of a fast and a slowly 
varying function of z. With the help of Eq. (|15|) one can 
write Eq. l(Ti|l as 



c 2 e iuz/c | y 



c az 



Let us now write the equation for the sources with the 
help of the curvilinear abscissa s, simply defined as s = 
vt. Here we assume that v =| v(t) | is a constant. We 
will use the general property 



*[/(*)] = £ 



I /'fo) I 



(17) 



where the zeros of f(x). Then, we will regard 

S(z — z'{t)) as S[f(t)], for any fixed value of z. The only 
zero of f(t) will be readily indicated with i(z), which is 
the value of t for which z — z'{t) = 0, while | f'(t(z)) \= 
v z (z(tj). Therefore, with the help of Eq. ltT7|) we can 
write 



5(z-z'(t)) 



V z (z) 



5(t t(z)) . 



(18) 



Then, using s(z) = vt(z), Eq. © and Eq. (gj give 



p(r ± ,z, t) = £ —S (r ± r' ± (z)) s( S -^-t) (19) 



and 



(r x ,z, t) = e —Mz)5 {r ± - r' ± (z)) 8 - t 



3 

so that 



p(rx,z,u) 



v z {z) 



(20) 



S (r± - r' ± (z)) e^ s{z)/v (21) 



and 



j{r±,z,uj) 



v(z)S (r_L - r' ± (z)) e 



ius(z) j v 



By substitution in Eq. l|lCj|l we obtain 



(22) 



2iu> d 
c dz 



Aire 

v z (z) 



cxp < iuj 



S[Z) 
V 



v(z)6 (n - r'Jz)) - VS ( Pl - r' ± (z)) 



(23) 



Eq. lIl'MI) is still fully general and may be solved in any 
fixed reference system (x,y,z) of choice with the help 
of an appropriate Green's function. In general, E will 
vary, in z', on a characteristic length which is related 
with exp {ico {s/v — z/c)}, which enters in the Green's 
function solution of Eq. 123(1 as a factor in the integrand. 
As it will be clearer after Paraerraph lll CI as we integrate 
along z' , the factor u>(s(z')/v — z'/c) grows larger and 
larger leading, eventually, to a highly oscillatory behavior 
of the integrand which does not contribute anymore to 
the final integration result. The value of z' for which 
ui(s(z')/v — z'/c) ~ 1 can be considered as a measure of 
the value of z' for which the integrand starts to display 
such oscillatory behavior and it is naturally defined as the 
radiation formation length L f of the system at frequency 
w. It is easy to see by inspection that if v is sensibly 
smaller than c (but still of order c), i.e. v ~ c but I/7 2 ~ 
1, then Lt ~ A. On the contrary, when v is very close 
to c, i.e. I/7 2 <C 1, the terms ujs{z')/v and —luz'/c tend 
to compensate so that Lf ^S> A. The exact expression for 
Lf depends on the particular magnetic system, i.e. on 
s(z'), and on the frequency of interest lu. 

Let us consider the assumption I/7 2 <§; 1. We will hold 
this assumption fulfilled throughout our paper. In gen- 
eral, introduction of a small (or large) parameter in any 
theory brings simplifications. In particular the ultra rela- 
tivistic approximation has two consequences: first, as just 
remarked, the radiation formation length is much larger 
than A and, second, as it will be discussed in Paragraph 
III CI radiation is emitted in a narrow cone of opening 
angle much smaller than unity. Accounting for these fea- 
tures of the radiation, and considering an electron mov- 
ing on a given trajectory, there will always be, in prac- 
tical cases of interest, some privileged reference system 



FIG. 1: Geometry and reference system for Synchrotron Ra- 
diation problems in paraxial approximation. 



FIG. 2: Standard reference frame for Synchrotron Radiation 
computations. 



siderably simplifies Eq. (|2"5|) as 



(x, y, z) in which Eq. (|23[l is simplified in a paraxial form, 
due to the ultra relativistic assumption, for some set of 
observer points. 

To apply such paraxial approximation we should first 
assume ultra relativistic motion and identify a main ob- 
server P. Then we can define a frame with the z axis 
oriented as the tangent from some trajectory point T, to 
P. As it will be discussed in Paragraph III CI using this 
frame we can describe radiation at P from all points of 
the trajectory with velocity vector forming with the z 
axis an angle much smaller than unity. In other words, 
for all these points 7^ 3> 1: this condition is necessary 
and sufficient for the applicability of the paraxial approx- 
imation. This is not a restriction, in practice, because it 
still keeps open the possibility of angles much larger or 
much smaller than 1/7 (bending magnet or undulator 
case) . The situation is described with a particular exam- 
ple for the case of a circular motion, in Fig. Note that 
the reference frame used to describe radiation at P can 
be also used for observers in the neighborhood of P, pro- 
vided that these form an angle much smaller than unity 
with respect to the velocity at T. Since we presented, 
as particular example, the case of a circular motion, it 
is worth to underline the difference between the frame 
depicted in Fig. ^and the standard frame used for Syn- 
chrotron Radiation computations in Fig. [21 using the 
frame in Fig. Elwe will never be able to describe the field 
at to observer locations displaced along the x axis. This 
is required, for instance, if one needs to calculate the au- 
tocorrelation function of the field at those two points. 
On the contrary, such a description is allowed with the 
choice of a frame like in Fig. ^ 

Since the radiation formation length is much longer 
than the wavelength, E does not vary much along z on 
the scale of A, that is | d z E XyV |<§; ui/c \ E XyV |. Therefore, 
the second order derivative with respect to z in the V 2 
operator on the left hand side of Eq. (|23[1 is negligible 
with respect to the first order derivative. This means 
that we can apply a paraxial approximation which con- 



l -^v ± (z)6 (r_L - r' ± (z)) - V ±6 (r x - r' ± (z)) 



2ioj d 
c dz 



Aire 



■ exp < ilu 



s(z) 



, (24) 



where we consider transverse components of E only and 
we substituted v z (z) with c, having used the fact that 
V7z ^ 1- Eq. 1|24[1 is Maxwell's equation in paraxial ap- 
proximation. Note that this approximation transformed 
Eq. ll'MI) which is an elliptic partial differential equation, 
into Eq. (|24l) . which is of parabolic type. 

The Green's function for Eq. (|24[1 . namely the solu- 
tion corresponding to the unit point source, satisfies the 
equation: 



V ± 2 + — G(z a - z; r ±a - r ± ) 

c oz ) 

= 5(r± - r ± )6(z - z) , (25) 

and, in an unbounded region, can be written explicitly 

as 



G(z Q ~ z';r ±c 



x exp < ilu 



r± 



4tt(z - z') 

|2' 



2c(z - z<) 



(26) 



As usual we will denote source coordinates with 'primes 
and observer coordinates with the index o. With the aid 
of the Green's function G, the solution of Eq. (|24|l can 
be represented as 



Ej_(z ,r± Ol uj) 



dz'- 



dr' ± 



^v ± (z')S (r' ± - r' ± (z>)) - V' ± S (r' ± - r' ± (z')) 
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x exp < iuj 



r± 



2c(z - z') 



s(z>) 



(27) 



where represents the gradient operator with respect 
to the source point. The integration over transverse co- 
ordinates can be carried out leading to the final result: 



~ . . live 
Ej_(z ,r ±0 ,L}) = T 



dz'- 



v x (z') 



Xq x (^z ^ 



X 



Vy(z') Vo ~ y'(z') 



Zq Z J \ C Zq Z 

where the total phase <!>t is given by 



(28) 



$7 



s(z') 



{x -x'{z')f + { Vo -y'(z')f 



2c(z - z') 



(29) 

Eq. (|28|l can be used in all generality to characterize radi- 
ation from an electron moving on any trajectory as long 
as the ultra relativistic approximation is satisfied, and 
a reference system suitable for paraxial approximation 
exists, which is always the case in situation of practical 
interest. Note that both x and y polarizations terms in 
Eq. (|28(1 are a sum of two parts which can be traced 
back to current and charge densities: in fact, the first 
part, proportional to the particle velocity in the x or 
y direction, follows from the transverse current density 
j±. The second instead, follows from the gradient of the 
charge density Vp. Our result makes a consistent use of 
paraxial approximation and of harmonic analysis which 
brings simplicity and power to the method. In the fol- 
lowing Sections we will show how Eq. (|28|) can be used 
to characterize bending magnet radiation and undulator 
radiation. 



xexp{z^| r a -r' |} . (31) 



Use of explicit expressions for p, Eq. 
(|22|l . leads straightforwardly to 



v z {z>) 



E(r ,u) = 
(3 — n ic 



icue 
c 



(|2T|) . and j, Eq. 



dz' 



n 



x exp < iuj 



r'{z')\ 

s(z') 



to |r -r'(z')| 2 
\r -r'(z>)\ 



(32) 



Finally remembering z' = v z t' we obtain Eq. (0). This 
result simply confirms that Eq. J2J is derived under most 
general conditions; the derivation that we proposed from 
direct use of fields is much less involved, from a logical 
viewpoint, with respect to the original, which makes use 
of potentials (see |2j||3(|) and it may be interesting for 
educational purposes. 

It is interesting to note that the only approximation 
applied to Eq. i|14|) in order to obtain, finally, Eq. 1281) 
was simply the paraxial approximation. We can easily 
see that if we apply paraxial approximation to Eq. @ 
we get back Eq. (|28|) . Let us show how this is possible. 

The choice of a; depends on our interest. Then, our ul- 
tra relativistic approximation implies a formation length 
Lf ^> A. Moreover in practical cases, \r a — r'\ will always 
be at least of order Lf, so that |r Q — r'\ 3> c/uj and, with 
an accuracy A/(27rL/), Eq. (0 can be simplified as 



E(r a ,u) = 

c 



dt'- 



(3-n 



\r a - r'(t> 



x exp \iu:[t'+ |r ° r/(f/)l ) \ (33) 



or, using again z' — v z t' ', v z ~ c and also s = vt' 



C. Discussion 

Following the previous derivation very natural ques- 
tions arise which regard the relation between our expres- 
sion, Eq. J5SJ), with Eq. Q and Eq. J5J) as well as the 
applicability region and the accuracy of our Eq. J2HJ- 

In order to investigate these subjects we go back to 
the most general Eq. (|14f) and we note that it can be 
solved by direct application of the Green's function for 
the Helmholtz equation: 



G{r a 



4tt 



* ex P {^|r o -r.'|}.(30) 



Integrating by parts the term in Vp we have 



E 



dr' 



C TV 



pn 



pn 



E(r ,u) = — 



dz 



, P-n 



x exp < iuj 



s(z') , |r -r'(0 



(34) 



Remembering that E± = E± exp{iu;z /c} we write the 
transverse field components as 



S i \ i^e [°° , (3± - n± 



x exp yiuj 
\r -r'(z') 



s(z') 



V c 

Zn - Z' 



(35) 



In any situation, it is mathematically correct to expand 
\r a — r'\ as an infinite sum of terms 



|r -r'| ~ {z -z') + 



\r_Lc 



2{z - z>) 



(36) 
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Truncation of the series in Eq. (|36|l and subsequent sim- 
plification of the phase in the integrand of Eq. 135|) is a 
delicate business. 

The first term of the expansion Eq. (|36[) naturally 
cancels the term — iuj(z ~ z')/c in the phase of Eq. 
iPH so that one is left with phase contributions due to 
higher order terms from the expansion Eq. 1|36|) and with 
iu>(s(z')/v — z'/c). This last contribution depends on 
s(z') and can be specified only when the magnetic sys- 
tem is specified. As has already been said in Paragraph 
III Bl as we integrate along z', uj(s(z')/v — z'/c) grows 
larger and larger leading, eventually, to a highly oscilla- 
tory behavior of the integrand which does not contribute 
anymore to the final integration result. The value of z' 
for which lo(s(z')/v — z'/c) ~ 1 can be considered as a 
measure of the value of z' for which the integrand starts 
to display such oscillatory behavior and it is naturally de- 
fined as the radiation formation length Lf of the system 
at frequency u>. 

Another length dictated by the physics of the problem 
is the natural size of the system. We will refer to it as the 
characteristic length of the system and we will indicate 
it with L c h', for instance, in the case of a circular motion, 
L c h is simply the circle radius. 

Relation between L c h and Lf depends on the system. 
In the ultra relativistic approximation we can say that 
L c h can be either much larger or comparable to Lf, but 
in any case never smaller, since integration of the Green's 
function is performed along a path length comparable 
with L c h- Moreover, in the ultra relativistic approxima- 
tion the two terms ujs(z')/v and - u>z'/c nearly compen- 
sate leading to a formation length much longer than A so 
that Lf ^S> A, and therefore also L c h ^> A. It should be 
stressed that, although the previous properties are very 
general, the relation between Lf, L c h and A depends on 
the particular physical situation under study and can be 
better specified only in relation with that situation: for 
instance, in the case of ultra relativistic motion on a cir- 
cular trajectory Lf — [\R 2 / (2n)] 1 ^ 3 . Also, note that 
if the ultra relativistic approximation cannot be applied 
anymore, then even the general scaling laws between L c h, 
Lf and A change. For instance, when if v ~ c but 7 ~ 1, 
one has L c h ~ Lf ~ A. 

The expansion Eq. (|36|) makes sense only if uj(s(z') fv— 
z'/c) is smaller or comparable to unity since the two 
terms us(z')/v and —luz'/c can be grouped together in 
a useful way. Once a certain wavelength of interest is 
fixed, the condition lu(s(z')/v — z'/c) ~ 1 determines the 
formation length Lf, which can be used in the second 
term of Eq. H36|) . By imposing that also this second 
term is not larger than unity one obtains a condition on 
the observation points of interest: 



\r_Lo 



2(z - z' f 



< 



2ttL, 



< 1 



(37) 



where we assumed z Q ~ z' > Lf. The third term in the ex- 
pansion can be then neglected, together with second and 



next terms in the expansion of the denominator of the 
integrand in Eq. Q35|). with an accuracy A/(27rLy). Note 
that, whatever the magnitude of r' ± across the beamline, 
it follows from Eq. I|37ll that the square of the opening 
angle of radiation r\ /z 2 , is much smaller than unity. 
This justify our attention to the transverse components 
of E only. 

If we substitute the first two terms of Eq. H36[) in the 
phase of the integrand of Eq. (|35[1 . and the first term in 
its denominator we find 



Ej_(r ,u) = -- 



s z 
x exp 



00 dz '^- n - 



\r± c 



Zo ■ 
12 



2c(z 



(38) 



that is just Eq. (|2g)> . 

We demonstrated that application of paraxial approx- 
imation to Eq. J2J) gives us back Eq. This fact 
is, however, quite obvious. In our derivation of Eq. i|28|) 
we applied paraxial approximation to Maxwell equations 
from the very beginning, while here we simply apply the 
same approximation after solving Maxwell equations in 
their more generic form. For consistency reasons we had 
to obtain the same result. 

The previous discussion tells that our approach is a 
specialization of the more general approach described by 
Eq. The novelty in this, stems from the fact that 

Synchrotron Radiation, by definition, applies to radia- 
tion by ultrarelativistic particles in magnetic structures; 
in Paragraph !!! Bl we demonstrated that paraxial approx- 
imation can always be applied to Maxwell equations in 
the case of ultrarelativistic particles. Eq. Q can always 
be reduced to Eq. I|28() in the treatment of Synchrotron 
Radiation. In other words, there is no point in starting 
from Eq. J2J) when treating Synchrotron Radiation prob- 
lems: one may start, directly, with the simpler Eq. i|28|) 
without loss of generality. 

Compared to Eq. J5J), Eq. (|28|l is easier to apply for 
analytical investigations, whose importance we already 
stressed in the introduction, because the paraxial ap- 
proximation has been applied from the very beginning, so 
that the approximations made, their region of applicabil- 
ity and their accuracy could be discussed independently 
of the magnetic system selected. 

The relation between our approach and the method 
in Eq. Q simplifies then, to the relation between Eq. 
and Eq. It is not straightforward to show this 

relation because Eq. J2J) does not coincide with Eq. 
once the term in \r a — r'\~ 2 is dropped. However, it 
is straightforward to show that Eq. J2Jl is completely 
equivalent to the expression for the Fourier transform of 
the Lienard-Wiechert fields including the velocity field 
part. This can be shown in full generality, regardless of 
the paraxial approximation. Eq. @ can be written as 



E(r ,u) = -e 



dt'- 



n 



r'(t'W 



x exp < iuj ( t' + 



dt' 



\ro-r'(t')\ 

c 

(3-n 



- exp < iu I t' + 



(l-n-/3)\r -r'(t> 
r -r'(t')\ 



dt' - l \ r 

where we have used relation 



(39) 



(4 



Eq. l(3"$l) can be integrated by parts. When the ed 
terms can be dropped (we will discuss this assumpti 
below) one obtains 



n 



K - r'(t> 



x exp < iuj 1 1' + 



dt' 



dt' 



(3-n 



(l-n- (3)\r -r'(t' 

. /, \r -r'{t')\\} 
x exp 1 iu it' + — — ] 



(41) 



With the help of Eq. igUJ) and 



dn 



dt' \r -r'(t')\ 



-—[-f3 + n{n-(3)} 



(42) 



one obtains, from Eq. 141fl the Fourier transform of the 
Lienard-Wiechert fields: 



E(r ,u) = -e dt' 



n (3 



7 2 (l-n-/3) 2 |r -r'| 2 



x exp \iuj\t' + ] — ' 



dt' 



1 n x [(n — /3) x p] 



(1 -n-/3) 2 



r ~ r'\ 

\exi)<{/V[t , + |r °~ r/(f,)l ) } . (43) 



The only assumption made going from Eq. I|39ll to Eq. 
(|43|l is that the edge term in the integration by parts is 
simply zero. This assumption can be justified by means 
of physical arguments in the most general situation ac- 
counting for the fact that the integral in dt' has to be 
performed over the entire history of the particle and that 
at t' — —oo and t' = +oo the electron does not con- 
tribute to the field anymore. It is obvious that the same 
line of reasoning can be followed starting from Eq. 143(1 
and going back to Eq. (|3^|) : in general, edge terms can 
be dropped. 

The previous statement is very general and per se triv- 
ial but it should be interpreted from a physical viewpoint 
depending on the situation. 



Screen 



B 




FIG. 3: A system constituted by two bending magnets con- 
nected by a straight section of length L. Radiation is collected 
at a distance z from the downstream magnet. 



For instance, consider the case of a system like the 
one sketched in Fig. 01 (which will be treated extensively 
in Paragraph IIII C|) . Fig. [3] shows two bending magnets 
separated by a straight line. Before and after the straight 
lines are two semi-infinite straight sections. A question 
arises about what we may take as t' = — oo and t' = +oo 
in this particular situation. Intuitively the magnets act 
like switches: the first magnet switches radiation on, the 
second switches it off. The statement that edge terms can 
be dropped, combined with the paraxial approximation 
tell that we can take t' = — oo and t' = +oo as the particle 
is well inside the first and the second bend, respectively, 
and neglect other parts of the trajectory. 

Note that an expression alternative to Eq. @ for the 
Fourier transform of the electric fields can be found in 
[3l| . After starting with Eq. 0, the authors of H3 
organized integration by part in a different way compared 
with what has been done in Eq. (|41l) . First they found 
that 



n x [(« - f3) x $] I d 

\r -r'\{l-n- P) 2 ~ \r a - r'\ ~dt' 



n x (n x p) 
(l-n-/3) 

n(n ■ p) + n(n ■ p) - n(n ■ pf - P{h ■ p) 



\r -r'\{l-n-pf 



(44) 



Note that Eq. 1(44(1 accounts for the fact that n = (r — 
r'(t'))/\r — r'(t')\ is not a constant in time. Second, 
using Eq. 1(44(1 in the integration by parts they got: 



E{r ,uj) = -- 



dt' 



n x (n x P) ic P — n — 2n(n ■ P) 
'\r -r'{t')\ + ~ 



x exp < iu> ( t' + 



|r - r'(t')l 2 
\ro-r'(t')\ 



(45) 



where, similarly as before, the edge terms have been 
dropped. Eq. 1(45(1 includes an integrand which is com- 
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pletely different with respect to Eq. (J5J). This is no mis- 
take. Both Eq. an d l|45l) are correct: using integra- 
tion by parts in different ways simply gives different in- 
tegrands which anyway, after integration, yield the same 
value for E(r a ,u>). It is interesting to note that, both in 
Eq. (J2J and Eq. (|45|l , terms in the first and second pow- 
ers of \r a — r'(t')| _1 cannot be interpreted as is usually 
done in time domain, as acceleration and velocity terms 
respectively nor they constitute, mathematically, Fourier 
pairs with E{r ,oj). In fact, a different organization of 
the integration by parts leads to different results: such 
terms have, here, no separate physical meaning. We have 
just shown that this situation is obtained starting from 
the Lienard- Wiechert field, applying integration by parts 
in two different ways and dropping the edge terms. 

Going back to the fact that edge terms can be dropped, 
we can make now some remark which is particularly in- 
teresting from a methodological viewpoint. Edge terms 
are important only when calculation of radiation proper- 
ties from a given system is performed over a part of it, 
that is integration is not taken from t' = — oo to t' = +00 
but only on a part of the trajectory arbitrarily chosen. 
However, summing up all non-negligible contributions 
is equivalent to integrate the system from t' = —00 to 
t' = +00. From this viewpoint edge radiation can be 
calculated without accounting for edge terms, which are 
artificial. As we will show in Paragraph IIII CI the low- 
frequency component of the field generated by a particle 
moving as in Fig. [3] usually referred to as edge radiation, 
actually arises from the straight line between the bends 
when one integrates the system over all the particle tra- 
jectory. 

As said before, on the one hand the term in \r a — r'\~ 2 
in Eq. @ can be dropped in paraxial approximation but, 
on the other hand, Eq. (|2J) does not coincide with Eq. 
Q once the velocity term is dropped. In other words, 
the term in |r Q — r'(t')| _1 of Eq. I|43|) is not the term in 
|r Q — r'(t')| _1 of Eq. @. Even in paraxial approxima- 
tion, the velocity term in Eq. (|43|l , in general, cannot be 
dropped. In fact, in order to do so, it is required that the 
ratio between the modulus of the velocity field and the 
modulus of the acceleration field is much smaller than 
unity: 



n- (3 



l 2 \r a 



r'\\nx [{n - P) x P]\ 
Then, in paraxial approximation we have 



< 1 



(46) 



nx[(n-f3)x$} = -nx[/3x(n-/3)} ~ [$-(n-j3)]n (47) 
Substitution in condition 1)46(1 gives 



7 2 |r G - r'\ 



n- f3 
n-f3 



P »1 



(48) 



Condition (|48|l depends on the system under investiga- 
tion, but in general is not automatically satisfied. For 



instance, in the case of circular motion /3/c 
being the circle radius. Then we have 



7>o 



P 



P 



R 



1/R, R 



(49) 



where 9 is the angle between n and 0. For instance, it 
is well known that the radiation formation length at the 
critical wavelength A c ~ -R/7 3 is simply Lf ~ R/j and 
that 9 ~ I/7; then, in the case r ~ L/, we have 



l 2 r o o 
R 



(50) 



As a result, when r ~ Lf and A ~ A c condition l|48|) 
is not satisfied, although the paraxial approximation is 
enforced. This counterexample shows that paraxial ap- 
proximation, alone, is not sufficient to guarantee that the 
velocity field in Eq. (|43|l can be dropped. 

The region of parameter space for which the velocity 
field can be neglected is usually referred to as the far field 
zone. In order to be in the far field zone the paraxial 
approximation is necessarily enforced. 

We take the condition that n is constant as a definition 
of the far field zone. When this is the case, the classical 
result (see for instance Q) 



n x [(n - P) x /3] _ 1 d 
(1 -n-p)\r -r'\ ~ T db> 



n x (n x P) 



1-n-P 



(51) 



can be used to perform integration by parts. Starting 
from Eq. (|43|l and assuming, as explained before, that 
the edge terms can be dropped, we arrive at the widely 
used result presented in textbooks (see again Q) 



E(r a ,u) 



zoje 



dt' [n x (n x p)\ 



x exp { iuj [ t' - -n ■ r' ) \ . (52) 



Now, with ultra relativistic accuracy 1 — n ■ P one has 
nx(nx/3) ~ P—n: using this relation we see straightfor- 
wardly that, in the far field zone, Eq. I|52|) is equivalent 
to Eq. H33fl . To conclude, in this discussion we have seen 
that paraxial approximation can be applied in the cal- 
culation of Synchrotron Radiation properties for any set 
of parameters. Moreover, in general, we should account 
for both velocity and acceleration term in the Lienard- 
Wiechert expression: however, in the case n can be taken 
as a constant (far zone approximation), the velocity term 
can always be dropped. 
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III. BENDING MAGNET RADIATION 
A. Circular motion 

Consider a particle moving along a circular trajectory 
and an observer as sketched in Fig. [I] The radiation from 
a particle moving with velocity v along a circle of radius R 
is observed from a point P. R is the only characteristic 
length of the system so that we can impose naturally 
that L c h is R. When 9 y and 9 X , defined in Fig. [I] are 
much smaller than unity the reference system [x, y, z) 
indicated in Fig. [J] is a consistent with the conditions 
for paraxial treatment explained in the previous Section 
and Eq. 128|) can be used to characterize the radiation 
at P. The motion along the curvilinear abscissa s can be 
described as 



r' ± (s) = -R(l- cos(s/R)) x 



(53) 



order in z' in the expression for s(z') in the phase, since 
the term z' /(2j 2 c) includes the small parameter I/7 2 . 

Note that, although integration is performed from —00 
to 00, the only part of the trajectory contributing to 
the integral is of order of the radiation formation length 
Lf = [XR 2 /(2t:)} 1 / 3 . At the critical wavelength i?/ 7 3 
that is simply ~ R/j- Physically, this is included in our 
paraxial approximation. Mathematically, it is reflected 
in the fact that <J>t in Eq. (|58|l exhibits more and more 
rapid oscillations as z' becomes larger and larger due to 
non linear terms: in particular note that, as z' ~ L/, the 
term in z' 3 in Eq. (|58l) is of order unity. If z a > Lf we 
can expand all expressions in (z a — z') -1 around z a . We 
will be considering the so-called far field radiation limit 
so that, expanding (z — z') _1 we will retain up to the 
third order in z' in the expression for the phase and up to 
the first order in z' in the rest of the integrand. We intro- 
duce angles 9 X — x /z and 9 y — y Q j z a as in Fig. [I] (with 
the restriction 9 X « 1 and 9 y <C 1). Then, accounting 
for v ~ c, Eq. """7|) and Eq. "Sl"ll read respectively 



and 



z'(s) = i?sin(s/r) 



(54) 



where s — vt 1 and the expansion in Eq. (|53() is justified, 
once again, in the framework of the paraxial approxima- 
tion. 

Since the integral in Eq. (|28|l is performed along z' 
we should invert z(s) in Eq. (|54|l and find the explicit 
dependence s(z'): 



s(z') = R arcsm(z'/r) ~ z' + 



6R 2 



(55) 



so that 
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r\ (z 1 ) = -—x 
x 2R 



(56) 



Substituting Eq. (55jl. Eq.JSU and Eq. (JS3J) in Eq. 
(12*51 we can write 



dz'e 1 ^ 



z' + R6 



R 



X A 

X 



(59) 



and 



<f>r = d> 



n + 



Zo + 



1 



q2 1 02 



2c "I ' V 2 7 2 c 2c 

I 6-r. \ n I 1 



x >z' 2 



,'3 



(60) 



\2Rc) ~ ' \6R 2 c, 
One can easily reorganize the terms in Eq. I|6U|I to obtain 



$ T = U! 



n + e 2 v \ 

2c Z ° 2c V7 2 ' 3 " 



ROx (± 
2c ^ 7 2 



-2+ d l 



(z' + R9 x ) (z' + R9 x y 



(61) 



2c 6R 2 c 
Finally, redefinition of z' as z' + R9 X gives the final result 



~ , . iLoe I ~~ , , e 
E±(z ,r± ,u) = / dz 



vz' + x + z' 2 /(2R) 
Rc z Q — z' 

where $t is given by 



— 00 



X 



Do 



7 V 



2 7 2 c 6R 2 v 



2c(z - z') 



(57) 



z' 3 [x + z' 2 /(2i?)] 2 + y 2 ^ 



At this stage our expression is still very general and valid 
for any observation distance z Q . We keep up to the third 



E=—< 

c A z Q 

x exp \^iui 
where 

and 



dz' ( — x + 9 y y 



2 7 2 c 



§o = — 



jR9. r ( 1 



2c V7 



6R 2 c 



*s = ^r(el + e 2 y ) 



. (62) 



(63) 



(64) 



Note that the linear term in (8 2 + 9y)z' in the phase of 
Eq. 159fl defines the maximal angle of interest. In fact, 
when z' ~ Lf, this term is of order unity if 9 2 + 9 2 ~ 
[A/(27ri?)] 2/3 ; this is simply the ratio of L 2 f to L 2 ch . As the 
observation angle grows more the linear term becomes 
larger and larger resulting in an oscillatory behavior of 
the entire integral in Eq. Q59JI an d thus giving no net 
contribution to the field. A similar mechanism has been 
discussed in relation with the radiation formation length. 
Now that we have found an upper limit to the observation 
angle, Lf / L c h, it is interesting to discuss the accuracy of 
our expression. We have 



axj> (rad) - x polarization component 



r±. c 



2(z - z'Y 



J2 



r 4 



T 2 



Rz AR 2 z 2 R 2 



< 1 . 

(65) 

In other words, the accuracy of our calculations is given 
by the square of the maximal observation angle. 

To compare Eq. I|62|) with results in literature (for 
instance |a,0) we should remember that the latter are 
obtained using a reference system like the one in Fig. 
Therefore, in order to perform a comparison, we should 
let 6 X — 0, so that $ = 0. Then, since results are often 
obtained in terms of integration along the retarded time 
t' we should use z' ~ vt' . Finally we obtain 



~ itoe iS) f°° , , / ct' 

CZo J-oo \ R 



x exp < ilu 



(66) 



Every Synchrotron Radiation textbook (see, for instance 
shows that Eq. (|6^|) can be written as 



CZ 07 J C 



#2/3 - i 



[1 + 7*02)1/2 



(67) 



Here K 2 /3 and K 1 / 3 are the modified Bessel functions of 
second kind of fractional order 2/3 and 1/3 respectively, 
while 



Ruj 
37 3 c 



(1 + 7^) 



(68) 



Note that, depending on the definition of Fourier trans- 
form Eq. I|12|) . textbook versions of Eq. I|66|) or Eq. 
(|67|l may differ by a factor 1 / ^/2tt and (or) by the oper- 
ation of complex conjugation; moreover textbooks often 
neglect the exact phase factor for the field, since they are 
usually more interested in the calculation of the intensity, 
which involves the square modulus of Eq. I|66|l . 

The result in Eq. (|67|l is far from being trivial. Bessel 
functions are real, therefore the x component in Eq. I|67|) 




-0.01 



FIG. 4: Phase correction for the horizontal polarization com- 
ponent of the field; case discussed in [27j . 



is real, while the y component is purely imaginary. In 
other words, the x component of the integral in Eq. I|66|) 
is purely imaginary, while the y component is real. This 
is not immediately obvious, but it can be easily seen by 
inspection, accounting for the fact that the exponential 
function in the integrand can be written as exp(i • ) = 
cos(-) + isin(-); then, for parity reasons, the term in x, 
being odd, couples with the sine function, thus giving an 
imaginary result while the term in y, being even, couples 
with the cosine function, giving a real result. 

From the previous remarks it follows that the only 
non-trivial phase factor is specified by the surviving ex- 
ponential argument <& s — ujz o 2 / '(2c), which is usually 
neglected in literature. This simply represents, in our 
paraxial approximation, the phase difference between the 
point (0, y , z ) and the point (0, 0, z )- Physically, in the 
particular reference frame of Fig. [21 which we have cho- 
sen by setting 9 X = 0, the electric field is represented by 
a spherical wave propagating outwards from the origin of 
the coordinate system. 

It is interesting now to investigate the meaning of 
the phase <I> , which is nonzero for nonzero values of 
8 X . This term has no equivalent in the usual treatment 
of Synchrotron Radiation, because displacement along 
the x axis cannot be considered in the reference system 
sketched in Fig. [2] In the usual treatment, an x dis- 
placement is treated redefining the reference frame so 
that x — 0. This redefinition has simply the effect of 
shifting the phase of the field expressed in the old refer- 
ence frame of a quantity i$ . 

There is a considerable amount of literature dealing 
with the phase of the Fourier components of the elec- 
tric field vector in Synchrotron Radiation (see, for in- 
stance H3). In particular, [^J reports the results of a 
thorough simulation of single particle effects. It is obvi- 
ously found that the wavefront of a single particle is not 
spherical; a correspondent phase shift (with respect to 
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ax> +7c x(9)(rad) - y polarization component 




y (m^o.oo5 



x (m) 



-0.01 



FIG. 5: Phase correction for the vertical polarization com- 
ponent of the field; case discussed in [2^1. 



of the statistical properties of a light source. We have al- 
ready discussed the advantage of our method regarding 
calculations of the field autocorrelation function in con- 
nection with the choice of a fixed reference frame as in 
Fig. ^ Let us now discuss how it can be used to calculate 
E from a single particle with a given angular deflection 
with respect to the orbital plane of a nominal electron. 
Such an expression was first calculated, starting from the 
Lienard-Wiechert fields, in [32|. 

The meaning of horizontal and vertical deflection an- 
gles rj x and rj v is clear once we specify the particle velocity 



v(s) 
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- cos 



cos{ij v 
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- sin(?? y )y 
cos(r] y )z 



(69) 



so that the trajectory can be expressed as a function of 
the curvilinear abscissa s as 



the spherical wave) for the field of a single particle was 
computed and analytical estimations within 1% of the 
numerical result were presented which are in complete 
agreement with Eq. Q64)). For the sake of completeness 
and comparison we present, in Fig. ^ and Fig. [5] the 
phase shifts calculated for the horizontal and vertical po- 
larization components, respectively by means of our Eq. 
(|64|l for the same example considered in j23|: bending 
magnet emission from a 2.5 GeV particle, constant mag- 
netic field in bending magnet 1.56 T, photon energy 40 
eV and distance from tangential source point to optical 
component 5 m. Note that these parameters correspond 
to the far zone case where z a ^> R/'y. The phase for the 
?/-polarization component, in Fig. [5] is simply the phase 
for the x-polarization component (that is Fig. ^) added 
to nH(0 y ), where H(6 y ) = 1 for y > and H(6 y ) = 
for 6 y < 0: this last terms simply accounts for the fact 
that the y component of Eq. lj6U)l is odd. 



l x + Rcos \ — + rj x ) cos(r)y) 
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Here we have introduced, also, an arbitrary offset 
(l x ,l y ,0) in the trajectory. Using Eq. (|7U|l an approx- 
imated expression for s(z') can be found: 



so that 
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B. Circular motion with offset and deflection 

Electrons following a circular motion are, of course, an 
approximation. Electron beams have always some small 
angular spread and offset with respect to the nominal tra- 
jectory. Any beam with a small geometrical emittance 
can be thought, in agreement with paraxial treatment, 
as a composition of perfectly collimated beams with dif- 
ferent deflection angles with respect to the orbital plane 
of the nominal trajectory. This representation is use- 
ful when one is interested in calculating, for instance, 
the influence of angular spread and offset on the electric 
field intensity: one can simply compute the contribution 
for each collimated beam and sum up the results. This 
method allows to simplify calculations of more compli- 
cated quantities like the field autocorrelation function, 
which is of uttermost importance in the characterization 



r' ± (z') = + r] x z' + l x ^x + (riyz' + l v )y . (73) 

We will consider ^ C 1 and r\ y <C 1 in agreement with 
the fact that we deal with electron beams with small an- 
gular deflection and not with more general kind of plas- 
mas. 

Substituting, as in the previous Paragraph, Eq. I|71(l . 
Eq. (O and Eq. (J22J in Eq. we can write 
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Redefinition of z 1 as z' + R(0 X — rj x ) gives the result 



where $t is given by 
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At this stage, as in the previous Paragraph, our expres- 
sion is still very general and valid for any observation 
distance z a , and again as we did before, we have retained 
up to the third order in z' in the expression for s(z') in 
the phase, since the term z' /{2j 2 c) includes the small 
parameter I/7 2 and up to the first order in z' in the 
rest of the integrand. We will now consider the far field 
radiation limit which, as the reader remembers, allows 
expansion of all expressions in [z a — z')~ l around z : as 
before we will retain up to the third order in z' in the 
expression for the phase and up to the first order in z' 
in the rest of the integrand. From Eq. (|74*jl . Eq. {TSjl 
and Eq. I|73|) it is evident that the offsets l x and l y are 
always subtracted from x Q and y respectively: shifting 
the particle trajectory on the vertical plane is equivalent 
to a shift of the observer in the opposite direction. With 
this in mind, in analogy with Fig. ^ we introduce angles 
Ox = X — lx/zo and 8 y = 8 y — l y /z (with the restriction 
9 X <C 1 and 9 y <C 1). Then, accounting for T) x <C 1 and 
r] y <C 1, Eq. (f7^|l and Eq. l(75|l can be written down 
respectively, as follows: 
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Except for the phase term $ s , Eq. (|79|l can be obtained 
from Eq. H62[) simply by substituting 6 X with 9 X — 77^ 
and with 9 y — rj y : besides a common phase factor 
then, we can say that including a deflection angle has the 
same effect of shifting the observer position of the same 
angle. Still, we should remember that 9 X = 8 X — l x /z 
and 9 y = 9 y — l y /z a . 

In the limit for l x , v /z < 6 F , with 6 F = [A/(27ri?)] 1 / 3 , 
one can simplify further Eq. I|79(l . Note how introduction 
of small parameters allows increasing specialization and 
simplification of the theory: for instance, in this Para- 
graph we started with an expressions for the fields ob- 
tained in the limit for I/7 2 ^ 1, that could be cast into 
simpler form in the limit Lf/z <C 1 and now further 
specialized assuming l x . y /z <C Op leading to: 
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One can easily reorganize the terms in Eq. (|77|l to obtain 
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Eq. H79fl is an extremely useful tool, because it describes 
the radiation from an electron with offset and deflection 
as in an electron beam with finite emittance, including 
the correct phase factor for the field. Starting from Eq. 
(|79|l then, it is possible to calculate the field correlation 
function, and to provide a study of transverse coherence 
properties of the radiation from a given electron beam by 
means of analytical techniques. 
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As has already been remarked, coherence is a spe- 
cial and important property of Synchrotron Radiation 
sources. In the far zone, the spatial coherence at differ- 
ent observer angles 9 x , y and fll 2 ,, is characterized by the 
autocorrelation function (see [3y, |33| and [38| for appli- 
cations to Synchrotron Radiation science): 



r (e (1) 9 (2) i 



= (E 



where Ej is the i-th Fourier component of the electric field 
given in Eq. {75J), and brackets (...) denote an ensemble 
average with respect to electron parameters. Note that 
a given electron is correlated just with itself, that is why 
Eq. (|85|l only includes the electric field from a single elec- 
tron and not from two different particles. The ensemble 
average can be replaced by integration over the electron 
phase-space density. The intensity distribution can be 
obtained directly from Eq. (|85|l by letting 9x}y — 6x}y 

Of course a numerical code can always be developed, 
either starting from Eq. i|28fl or just from the Lienard- 
Wiechert fields, which calculates the field correlation 
function in a generic case, but such a code would not 
help in physical understanding of the situation. On the 
contrary, Eq. (|79|l includes all relevant information about 
an electron in a realistic beam (i.e. with offset and deflec- 
tion) and, being an analytically manageable expression, 
constitutes the first step towards the characterization of 
transverse coherence properties from bending magnet ra- 
diation. 



C. Edge radiation 

In this Paragraph we will stress the importance of the 
knowledge of the entire trajectory followed by the elec- 
tron as we study the effect of a change in longitudinal 
velocity due to the passage of an electron in a magnetic 
system. This results in collimated emission of radiation 
in the low photon energy range, a mechanism analogous 
to transition radiation, which is well-known in literature 
under the name of edge radiation [53, H^, 0] . 

We restrict ourselves to the system depicted in Fig. 
[3J which shows two bending magnets separated by a 
straight line of length L. Radiation is detected at a 
screen positioned at distance z from the downstream 
magnet. We will require that the bending magnets de- 
flect the electron trajectory of an angle much larger than 
[A/(27ri?)] 1 / 3 , A being the wavelength of interest: in this 
way the straight lines before the upstream bend and after 
the downstream bend do not contribute to the field de- 
tected at the screen position. Intuitively, the magnets act 
like switches: the first magnet switches the radiation on, 
the second switches it off. The trajectory can be there- 
fore split in three parts: the two bends and the straight 



section of length L between them. With the help of Eq. 
(|28[l we write the contribution from the straight line as 
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where we assumed z a 3> L. The previous assumption is 
not always verified in cases of practical interest. Here, 
however, we are only concerned with an example of ap- 
plication of our method. <&t in Eq. i|86|) is given by 
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The characteristic length related to the straight section 
is obviously L c h = L. In general, the formation length of 
radiation at wavelength A is given imposing ioz 1 j (2c7 2 ) ~ 
1, which gives Lf ~ 7 2 A. However, if we are interested in 
low frequencies such that ujL/(2j 2 c) <C 1, we can simply 
consider Lf ~ L and neglect the term in I/7 2 in Eq. 
(187(1 . Trivial calculations show that 
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where $ s has the usual meaning (compare, for instance, 
with Eq. (|63|l % ). The magnitude of the contribution from 
the straight section AB can be estimated noting that the 
sine function in Eq. drops rapidly as 9 2 -\-0 2 reach the 
value Ac/(loL) = 2\/(ttL). This means that a maximal 
observation angle of interest related to the straight line 
can be found in the transverse direction: 
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At this angle, the formation length for the straight line 
radiation at wavelength A is simply equal to the straight 
section length, L, as it can be seen from Eq. H87fl . Then, 
the magnitude of Eab is of order 
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On the other hand, the magnitude of the contributions 
from the bends can be estimated as 
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In Eq. I|91|) we have redefined the origin of our reference 
system, as it is clear from the integration limits; however, 
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here we are interested in the magnitude and not in the 
phase of the field, so that we can use Eq. (|91|l without 
further discussion on the correct phase. The linear term 
in z' in the exponential function in Eq. (|91|l is of order 
unity for angles 9 y such that 

-* 9 * 1 ■ 



(92) 



Substituting the formation length for the bends LfB ~ 
[Ai? 2 / (271-)] 1 / 3 we obtain a condition for the maximal ob- 
servation angle of interest related to the bending magnets 
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y \2ttR, 

The previous condition on # 2 found for the bends has to 
be compared with the condition on # 2 previously found 
for the straight line, (9 2 ~ X/(2wL). The ratio between 
the latter and the former is equal to the ratio of the 
formation length for the straight line to the formation 
length for the bends respectively that is L/LfB- As- 
suming L/Lfg 1, that is always verified in practice, 
we obtain that the maximal observation angle of interest 
related the straight line is much smaller than the maxi- 
mal observation angle of interest related to the bend. As 
a result the entire linear term in the phase of Eq. 191(1 
can be neglected thus leading to the following simplified 
expression for E^: 
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An estimate of Eq. $Q is 
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The ratio between the magnitude of the bend and the 
straight section contribution is then 
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Therefore, under the already accepted condition L ^S> 
Lf one can neglect the contribution from the bending 
magnet. The only remaining contribution is given by Eq. 
(|88p. which represents an expression for the field at the 
detector position. We plot the intensity corresponding to 
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FIG. 6: Intensity distribution of radiation from the setup in 
Fig. E]at A = 10/im and for L = 1 m. Radiation is collected 
at z = 10 m from the downstream magnet. A contour plot 
(upper figure) and a horizontal cut by the median plane (lower 
figure) are presented. 



the field in Eq. in Fig. HJfor the case A = 10 //m, 
L = 4m and z = 20 m. 

Eq. iJSSJ is the same reported in 0> 0- It is in- 
teresting to remark that in conventional treatments of 
Synchrotron Radiation, edge radiation has its origin in 
the edge term arising in the integration by parts of Eq. 
||TJ ; in other words it is found from the acceleration part 
of the Lienard-Wiechert field, which is present only in 
the magnets. On the other hand, from our method, edge 
radiation arises as the contribution from the straight sec- 
tion between the magnets. This seems paradoxical, but 
one has to remember that there is no physical meaning in 
the calculation of the field harmonic content over a part 
of the trajectory alone. It does not make any sense, for 
instance, to calculate the field intensity from the edge of 
a single bend, because there will be either a second bend, 
as in Fig. |3 other structures downstream of the second 
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FIG. 7: Radiation field from an electron moving along a circle 
in the time and in the frequency domain 



FIG. 8: Radiation field from an electron moving along an arc 
of a circle in the time and in the frequency domain 



bend which must be taken into account in the calculation 
of the total field: the sum of all these contributions will 
give interference terms when the intensity is calculated. 
Then, accepting the viewpoint that only the knowledge of 
the entire trajectory of the particle brings physical sense 
to field calculations, there is no real contradiction: in our 
method edge radiation appears from the straight section. 
In the usual approach, instead, it appears from the edge 
terms in the integration by part of the acceleration field. 
Yet these terms, alone, have no physical meaning. This 
solves our paradox. 



radiation from short bending magnet 
( L m « R/7 ) 



time variation of electric field 
at detector in far zone 
( electric field averaged over time / ) 
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D. Short magnet radiation 

Our method can be used to compute radiation charac- 
teristics from a short magnet, characterized by a bending 
angle 2ip -C 1/7. This device has the quite interesting 
characteristic that the Fourier transform of the electric 
field in the far field limit has non-zero value for lo — > 
and that the critical frequency does not depend on the 
magnet radius but only on the magnet length. These 
features can be explained in a simple way. For circular 
motion the far field from a single particle has the prop- 
erty: 

Edt = . (97) 

-00 

When the electron moves in arc of circles with different 
angular extensions 2ij) though, the time-average of the 
electric field is nonzero; then, its Fourier transform has a 
nonzero component for lo — > 0. 

The first studies on long wavelength asymptotic can be 
found in |33j . For a review on the subject it may be in- 
teresting to consult 0] ; for a discussion dealing also with 
properties of CSR (Coherent Synchrotron Radiation) see 
[34| . In Fig. Fig. |H1 and Fig. we plot the energy 
spectrum of the radiation in the case of an electron mov- 
ing on circular motion, dipole magnet and short bend, 



FIG. 9: Radiation field from an electron moving along a short 
bending magnet in the time and in the frequency domain 



respectively. All figures refer to the x polarization com- 
ponent. As Eq. (|97|l does not hold we find a nonzero 
component of the energy spectrum at lo — > 0. As the 
magnet extension becomes smaller than I/7 the critical 
frequency depends only on the magnet length L m = 2Rip 
according to lo c ~ 2cj 2 /L m . 

The short magnet is a device of particular interest also 
because its radiation characteristics are related to the 
methodological issue introduced in Section [J From the 
very beginning, our method relies on harmonic analysis 
of the electric field, so that the entire trajectory of an 
electron is considered known. This gives no particular 
problem in the case of a closed trajectory, as discussed 
in the previous Subsections, but a question arises, now 
that a short magnet is considered, about the trajectory 
of our particle outside the short magnet. In general, the 
field at the observation position depends strongly on the 
trajectory followed by the electron before and after the 
magnet. So one should clarify what is the meaning of 
the term short magnet radiation by itself, when it is 
not specified what follows and what precedes the bend. 

Computing the field from the magnet alone as it is usu- 
ally done has a meaning, but only in a particular sense: 
one can always break up the beamline in several parts, 
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calculate contributions separately and finally add them 
up, accounting for the proper relative phase, to get the 
total field at the observer position. Then the short mag- 
net radiation is simply a part of the total field, a contri- 
bution to be added up to something else, which depends, 
case by case, on what precedes and follows the magnet. 

In conclusion, when we talk about radiation from a 
short magnet, we mean the contribution to the total field 
at the observer position P due to the short magnet alone: 
this makes sense as long as we consider it only a part of 
the total field, calculated separately just for computa- 
tional convenience. 

With this in mind we calculate the short magnet field 
from Eq. (|59|) limiting the integration to the magnet 
extent and using 2tp <C 1/7- This introduces a further 
small parameter in our system which leads to extra sim- 
plification of already known formulas derived for the case 
of a bending magnet with arbitrary extension. From Eq. 
(|59"|) we obtain 
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where $ s is defined as in Eq. In Eq. (|9"S1) we 

intentionally separated phase factors exp{iuj9 x z' 2 /2Rc} 
and exp{iuiz' 3 / (6R 2 c)} because we intend to neglect the 
second and expand the first around unity. For our sys- 
tem the characteristic length is obviously L c h — Rip. In 
general, the formation length of the radiation is fixed 
by the first phase factor, Lf ~ 7 2 A, but we will be 
interested in studying the long wavelength region such 
that Lf ~ Rip, that is the maximum attainable value, 
since it is always Lf < L c h- Then, the first phase 
factor imposes, simultaneously, loRipj (2 7 2 c) < 1 and 
ui9 2 Rip/c < 1. If on the one hand, from the first of 
these conditions uo max ~ 2j 2 c/(Ri/j), from the second 
condition follows that 9 2 < l/(27 2 ) and from the second 
phase factor one has oj9 x (Rip) 2 / (2Rc) ~ yip <C 1. On 
the other hand if, from the first of these conditions, 
u) -C 27 2 c I (Rip) , then from the second condition follows 
that 9 2 ~ c/(luRi/j) ^> l/(2y 2 ) and from the second phase 
factor one has u9 x (Rip) 2 / (2Rc) ~ (jip)/(y9 x ) < 1- 
It follows from these considerations that we can per- 
form the expansion of exp{iui9 x z' 2 /2Rc} up to the sec- 
ond term with an accuracy ('yip) 2 or better and neglect 
expjicjz' 3 / (6R 2 c)} with the same accuracy for any choice 
of lo. Then we can write 



Performing calculations and dropping negligible terms we 
find 
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If we define a function X[-Rif>.Rif>] which assumes uni- 
tary value on the interval [—Rip, Rip] and zero value else- 
where, then X[-Ri>,Rip] simply represents the magnetic 
field shape at the magnet position, and X[-Rif>,Rip] is its 
Fourier transform with respect to a. 

Suppose now the magnetic field has a generic shape 
B(z'). Keeping the same accuracy as before, we can gen- 
eralize Eq. llOOfl substituting relations 
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and 
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where m e is the electron rest mass, in place of x' = 
-z' 2 /(2R) and v x /c = -z'/R directly in Eq. JUJ. 
Performing approximations as before and integrating by 
parts two time the term in x 1 and one time the term in 
v x we find that Eq. (|100fl can be easily generalized. 

Results in literature are often given in cylindrical co- 
ordinates ((f>, 9, z) where sin</> = y Q / \J x 2 + y 2 and 9 = 
9 2 + 9 2 . The same choice can be found, for instance, 0. 
Making use of this system we obtain: 
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where r e is the classical radius of the electron 
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and 



(106) 
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As has been said, the field from a short magnet makes 
sense only as a part of the total field from a given beam- 
line. In this Paragraph we have calculated, for such con- 
tribution, the expression in Eq. (|105fl which can also be 
found in [7J. Our aim, here, is not to review a well-known 
result, but to cast new light in its physical interpretation. 

In order to do so we need first a digression. Our Eq. 
(|28|l is based on the more general Eq. @. In Para- 
graph [^21 we have shown that Eq. (|2J) is equivalent to 
the Fourier transform of the Lienard-Wiechert fields Eq. 
(|43JI provided that the edge term in the integration by 
parts going from Eq. (|43|) to Eq. @ (or viceversa) can 
be dropped. On the one hand we showed that this can 
always be done since the integral in Eq. (|43|l and Eq. 
@ should be performed over the entire electron history, 
and only a finite part of the trajectory contributes, prac- 
tically, to the electric field at the observer position. On 
the other hand, one is free to break up the beamline in 
parts and sum up partial contributions to the total field 
using Eq. (|43|l followed by integration by parts for each 
segment. In this case one must retain edge terms to ob- 
tain the same result: as an example of this fact, using 
our method, we have shown in Paragraph IIII Cl that the 
edge radiation from the system depicted in Fig. [3] arises 
as the contribution from the straight section between the 
magnets, while conventional calculations indicate its ori- 
gins in the edge term from the integration by parts of Eq. 
(|4*3j) in the far field limit. 

Now coming back to the subject of this Paragraph, typ- 
ical derivations of Eq. ((TUB) (see USUI) start from 
Eq. Q and include the complete expression for the accel- 
eration field. Then, following conventional derivations, it 
is not clear wether edge terms play an important role in 
Eq. (|105|l or not. Yet, we were able to obtain Eq. I|105|) 
without starting with Eq. ft}: in fact we began with Eq. 
(|98J) which is Eq. I|59|) with different integration limits 
which is, in its turn, a reduction of Eq. (|28|) that, finally, 
is a simplification of Eq. (J2J. Since our result coincide 
with the one in literature it follows that edge radiation 
from the presence of the short magnet is completely ig- 
norable. 

Of course, one should still account for contributions 
from the other parts of the beamline: for instance, if the 
short magnet were installed in the middle of the straight 
section in Fig. |2|one should consider the edge radiation 
contribution as in Paragraph lHI CI 

However, according to our reasoning we can conclude 
that one can safely drop the contribution of edge terms 
due to the presence of the short magnet. On the contrary, 
beginning with Eq. Q as is done in usual treatments, 
one is not able to separate the two physical phenomena 
and, indeed, may be easily misunderstand and misinter- 
pret results concluding, erroneously, that short magnet 
radiation cannot be derived without edge contributions. 
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FIG. 10: Geometry for undulator radiation. 



Here we have seen that a critical study of the theoret- 
ical status of well-known formulas, can sometimes yield 
surprises. 



IV. RADIATION FROM INSERTION DEVICES 
A. Standard undulator 

For a review of up to date knowledge on undulator 
radiation one may be interested in consulting references 
|{| to [To) . An experimental characterization of radiation 
from insertion devices at third generation light sources 
can be found in [I^. Eq. (j2S|> can be used to derive the 
expression for E(r, ui) in the case of an undulator as well. 

Again, as before, we remark that, with the term undu- 
lator field we actually mean a part of the total field seen 
by an observer and that one should account, in general, 
for the entire motion of the particle. The situation is 
sketched in Fig. 1101 where a planar undulator consisting 
of N w periods is sketched. We will follow [5j in deriving 
well-known relations. For the electron transverse velocity 
we assume 



cK 

v±(z') = — sm(k w z')x . (108) 

7 

Here K is the deflection parameter, k w — 2tt/X Wi where 
X w is the undulator period, so that the undulator length 
is L w = N W X W . The undulator length L w is also, nat- 
urally, the characteristic length of the system, so that 
L c h = L w . The transverse position of the electron is 
therefore 



r\ iz') = —r- cos (k w z')x . (109) 

An expression for the curvilinear abscissa s as a function 
of z is given by 



8 K 2 
s(z) — — — z — — sin (2k w z') 



(110) 



8 av 8~f 2 k v 

where 8 av is the time-averaged velocity along the z di- 




FIG. 11: Radiation geometry for the undulator field. 



rection, that can be expressed as: 



FIG. 12: Undulator followed and preceded by bending mag- 
nets and straight sections. Information about the particle 
trajectory before and after the undulator is necessary to cal- 
culate the radiation spectrum. 



K 2 

'"' = 'I '"^ 



(111) 



We can now substitute Eq. (|109fl and Eq. I|1U8|) in our 
Eq. Such a substitution leads to a general expres- 

sion, valid for any observer distance z a . It is possible to 
obtain, similarly to many Synchrotron Radiation text- 
books, a simplified expression valid in the limit for large 
values of z Q . Since we are interested in the contribution 
of the undulator device to the total field at the observer 
position, we will integrate Eq. (|28|l only along the undu- 
lator. Then all terms in [z — z') _1 in the phase factor 
of Eq. I|28|) can be expanded around z Q . In the limit 
z Q 3> N W X W , we can usually retain first order terms in z': 
later on, in Paragraph II V CI we will discuss the applica- 
bility region of this approximation and we will see that 
there are particular regions of parameters which do not 
allow retention of first order terms alone. Dropping negli- 
gible terms and defining cj) and 9 according to Fig. ^] (or 
remembering the analogous definition given in Paragraph 
IIII D|) we obtain 



E ± (z 0l r± Ol uj) = — — / 

C Z a J-N w \ m /2 



dz 



K 



sin (k w z') + 9 cos (</>) I x + 9 sin (0) y 



(112) 



Here 



$ T = 



7 

K 2 W! . 



87 2 k w c 



sm(2k w z') 



$s , (113) 



where $ s = oj9 2 z /{2c), as usual, and u>\ is defined as 



2k w cj 2 V" ' 2 



K 2 2/1 . 
1 + — + 7 0' 



(114) 



Note that asw = wi, the formation length of our system 
is simply given by Lf = l/k w as it is seen immediately 
from Eq. (|113|) imposing the first term in phase to be 
of order unity. Yet it should be stressed here that the 
formation length Lf = l/k w is obtained in the most gen- 
eral case, without accounting for special conditions which 
are very often chosen for undulator operation. In fact, 
when a large number of undulator periods is selected an 
extra large-parameter is introduced in the system yield- 
ing simplifications within the theory; then, as it will be 

if the so called res- 



clear after reading Paragraph IIV Bl 
onance condition is met, the actual formation length of 
the system becomes Lf ~ L w . 

The choice of the integration limits in Eq. i|112|) ex- 
press explicitly the fact that the reference system has its 
origin in the center of the undulator as in Fig. ^| As said 
before, Eq. (|112fl gives only a partial contribution to the 
total field, which must be summed up with terms arising 
from structures preceding and following the undulator. 
In [5|, for instance, no particular attention is given to 
this problem: the calculation of undulator radiation dis- 
tribution is performed considering the acceleration term 
of the Lienard-Wiechert fields alone, integrating by parts, 
and dropping the edge term. 

It is indeed very useful, at this point, to comment some 
passage of 5]. At page 412 one can read: "No assump- 
tions on the magnetic field parameters have been made 
to derive the radiation spectrum in the form of equation 
(7.87) which we use to calculate the radiation spectrum 
from a wiggler mag net" . Eq. (7.87) in @ is the radiation 
spectrum from a particle moving on a generic trajectory 
calculated, indeed, starting from the Lienard-Wiechert 
field expression for the acceleration field, integrating by 
parts and neglecting the edge terms; it obviously appears 
in the form of an integration with limits from — oo to oo: 
this is just our Eq. I|52|l that was derived exactly as Eq. 
(7.87) in p|. Dropping the velocity field in the Lienard- 
Wiechert expression is done under the sufficient (but not 
necessary, in general!) condition that n can be considered 
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constant. This is regarded as the far field approximation. 
Then, neglecting of edge terms can be done under physi- 
cal assumptions discussed much earlier in Paragraph !!! CI 
Again, the integral in dt' has to be performed over the 
entire history of the particle; the physical meaning of 
t' = — oo and t' = +00 is that at these times, the elec- 
tron does not contribute to the field anymore. From this 
viewpoint, the geometry specified in Fig. 1101 is no more 
sufficient. One has to provide information about what 
precedes and what follows the undulator; an example is 
depicted in Fig. ^| In the case of the scheme in Fig. 
1121 for instance, t' = —00 and t' = +00 refer to moments 
when the particle is well inside the first and the second 
bend, respectively, and do not contribute anymore to the 
resulting radiation spectrum. In order to drop edge terms 
in the particular case of the trajectory in Fig. 1121 one has 
to consider contributions from the first bend, then from 
A to B along the straight line, from B to C inside the 
undulator, from C to D along the second straight line 
and, finally, from the second bending magnet. 

Instead of following this prescription, the author of Q 
substitutes characteristic quantities for a trajectory on 
an infinitely long undulator in Eq. (7.87), and then he 
simply sets the undulator strength parameter K to zero 
outside of a given temporal interval, which is simply the 
physical time that an electron takes to pass through the 
undulator. This is equivalent to consider two infinite 
straight lines before and after the end of the undulator. 
Then, although results are correct in some particular case 
(namely under resonant approximation, which we discuss 
in Paragraph llV Bfl the reader is left with the unsolvable 
task of understanding what is the physical meaning of 
two infinite straight lines and, especially, what is the con- 
tribution of these infinite parts of the trajectory to the 
total electric field at the observer position. Again, if one 
wishes to drop the edge terms, one should start consid- 
ering the entire trajectory of the electron, not a part of 
it alone. 

In the case depicted in Fig. ^| for instance, we will 
have different interfering contributions from the straight 
lines before and after the undulator and we would end up 
with both edge radiation and transition undulator radi- 
ation. Again, in the most general case, Eq. (|112f> refers 
to a part of the electric field seen by the observer and it 
should be added to other contributions corresponding to 
the entire history of the particle. Only in particular sit- 
uations like the one discussed in the following Paragraph 
II V B\ calculation of radiation properties from Eq. I|112|) 
alone makes sense: in Paragraph llV Bl we will discuss the 
resonant approximation which can be used in the case of 
a large number of undulator periods, N w 3> 1. Again, the 
presence of an extra large (or small) parameter yields to 
further simplifications of the theory. However it should 
be noted that, in practice, ranges from a few units, 
in the case of Far Infrared insertion devices to about one 
thousand, for X-ray FELs, so that the applicability of 
the resonant approximation strongly depends on the sit- 
uation under study. 



B. Resonant approximation 

In general, it does not make sense to calculate the in- 
tensity distribution from Eq. i|112|) alone, without extra 
interfering terms. Yet, we can find particular situations 
for which the contribution from Eq. (|112fl is dominant 
with respect to others. In this case, and only in this case, 
Eq. (|112f> . alone, is endowed with physical meaning. 

The method proposed in 0] to calculate the integral in 
Eq. (|112fl is very well-known and used by every textbook 
treating the details of undulator radiation. Again, we will 
follow notation in [j| during our derivation. The method 
consists in using the identity 
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where J p indicates the Bessel function of the first kind of 
order p. 

After introduction of 



Q = 



and 
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2K(3 avl 6cos(<P) 
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A{1 + K 2 /2 + j 2 8 2 ) 
one can express the exponential function in Eq. (|112ll as: 
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e ** T = E E Jm{u)J n {vy™l 2 e- m ^ z ' (118) 
where 



UJ LO UJ 

u= — S, v = — Q , and = n — 2m. (119) 

UJ\ UJ\ UJ\ 

Then, Eq. (jll2|) can be re-written in the more suggestive 
form: 



E±(z ,r ±0 , uj) = -^e l * s V" J m (u)J n (v) 
c z z r ' 



xe 



m : n= — 00 
N m X w /2 



inn/2 j ^ , I —iR^k m z' 

N w X w /2 



K 

— sin (k w z') + 8 cos (0) ) x + 8 sin (0) y 



Performing the integral in Eq. I|12l)[l one gets: 
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+ COS (0) 

+9 sin (0) 
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sin(7rAT lu i? ti; ) 
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As is well known, and can be seen by inspection, terms in 
Eq. l|121l) exhibit resonant character. Maxima of differ- 
ent terms are found as R u — 0, R u + l = and R u — 1 = 0. 
These correspond to particular frequencies multiples of 
the fundamental w\. Always following |{| we introduce 
the harmonic number k and we set ujk — ku)\. Then, 
after setting Auou — u) — uJh, the following expression can 
be obtained for the field around frequency Uk- 
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E±,k — 
K 



i<jjeA w N w e l<ia s'm(irN w Aujk /coi) 
c 2 z Q nN w AuJk/uJi 

t(J TO (tt)J fc _ 2 ,n-i(t;)e <w (*- 2m - :L 5/ a 



~J m {u)J k - 2m+ i(v)e^ k - 2m+1 ">' 2 ) 
+6cos( ( f>) J m (u)J fe - 2 „ l («)e l7r ( fe - 2m )/ 2 ] x 
+9sm(<P) J m (u)J k - 2m (v)e^ k - 2m ^ 2 y} . 



(122) 



The factor sin(-)/(-) represents the well-known resonant 
character of the device: our electron produces field 
peaked at frequencies io k only. u>k is a function of the 
angle 9. Once an observation angle 9 is fixed, a resonant 
frequency Uk (0) is also defined and the bandwidth of the 
radiation is determined by Alu^ (9) through the resonant 
factor sin(-)/(-). If we are interested in the angular width 
of the peak around the observation angle 9, we can intro- 
duce an angular displacement A9 with respect to 9. The 
resonant frequency at angle 9 + A9, that is Ldk{9 + Ad), 
has a different value with respect to 00^(9). When A9 
becomes large enough, uik (9 + AO) gets outside the band- 
width Auf-(0). This happens for an angular displacement 
A9 such that 



N w (uj k (0 + A9) - u) k {e))/u x {9) ~ 1 (123) 

which is the first zero of the sin(-)/(-) function in Eq. 
(|122ll . If we are interested in the angular width of the 
peak around 9 = in the case k = 1, that is for the 
fundamental harmonic, we should solve the equation 



N w ((Jx(Ae)-uJi(0))/cji(0) = 1 



(124) 



with respect to A9. Eq. 1|124|) should be taken only as 
a rough indication of the angular width of the radiation 
peak. In fact, taking the first zero of the sin(-)/(-) func- 
tion in Eq. (|122fl is just a convention: as we will see, 



only through similarity techniques we will be able to de- 
termine the natural angle which fits with the physical 
situation under study: 9 can be compared in a natural 
way only with that angle, and not with the solution of 
Eq. (|124f) . Yet, following textbooks and tradition we will 
call with 9 C the solution in A9 of Eq. Ijl24(l , and we will 
refer to it as a rough indication of the angular width of 
the radiation for the 1st harmonic. The cone with aper- 
ture 9 C is usually called central cone. It can be found 
that (see 0) 
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where 7 Z = y/a/1 + K 2 /2. From the previous discussion 
follows that Eq. I|122|) can be drastically simplified un- 
der the approximation N w 3> 1 and within the central 
cone 9 < 9 C . In particular when k = 1, that is for the 
fundamental harmonic, we have 



KojeX w N w e t ' t ' B sin(7r7V u; Aa; 1 /w 1 ) 
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Here we neglected all terms but the one proportional to 
N w . Also, Ajj{u) = Jo (w) — Ji (u), and we have used the 
fact that K < 7. 

Note that at this point we are not able to specify the 
accuracy of this approximation precisely: we can only say 
that, when 9 ~ 9 Cl our accuracy scales as 1/N W . This can 
be seen showing that the third term in x and the term 
in y in the series of Eq. (|122fl can be neglected with an 
error ~ l/N w . First we fix k = 1, that is the fundamental 
harmonic. Second we note, from Eq. (|119f) . that for any 
value of K and 9 much smaller than l/j z we have u of 
order unity while v is much smaller than unity, so that 
we can take it as a small parameter: then, the largest 
contributions from the sum in Eq. I|122|) will be for the 
smallest indexes of the Bessel functions J(v) because of 
the asymptotic behavior of J q (v) ~ v q . It follows that 
the third term in x and the term in y in the series of 
Eq. i|122D are of magnitude 9Jq(u)Ji{v) or 9J\{u)J-\{v) 
to be compared with the first two terms in x scaling as 
K/(2j)J Q (u)J {v) and K/(2j)J 1 {u)J (v). 

This situation refers to the case when trigonometric 
functions in (f> are of order unity, otherwise extra fac- 
tors cos 2 ((f)) or sin(</>) cos((f>) should be included as well. 
Taking the trigonometric factors of order unity gives an 
upper limit to the accuracy of our calculation and we will 
always do this. Whatever the value of u, neglecting the 
third term in x and the term in y in Eq. I|122|) can be 
done with an accuracy given by the ratio 
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Then, when 9 ~ 9 C we obtain i? ~ 1/N W . Yet, as has 
already been remarked, (9 C is not the natural angle to 
which 9 should be compared and as a result we cannot 
specify numerical factors in front of the 1 /N w . Later we 
will comment further on this point. 

Eq. 1)126(1 is indeed well-known, but usually no at- 
tention is paid to its deeper meaning. As we have said 
before, it does not make sense to calculate the inten- 
sity distribution from Eq. (|112|l alone, without extra 
interfering terms. Yet, we have told that, in particular 
situations, the contribution from Eq. I|112|l is dominant 
with respect to others and that in this case Eq. (|112|l . 
alone, is endowed with physical meaning. On the other 
hand, we have seen that Eq. (|112|l simplify to Eq. I|126|) 
if k = 1, N w ^S> 1 and within the central cone 9 < 9 C . 
Again, when 9 ~ 9 C such simplification is valid within an 
accuracy scaling with 1/N W . By inspection of Eq. Ijl26(l 
we see that, due to resonance, the only surviving term 
scales with N w . Any extra term added to Eq. I|126[l due 
to non resonant devices like bending magnets or straight 
lines will be simply negligible with respect the field in Eq. 
(|126[) with an accuracy scaling with 1 /N w . Note that Eq. 
(|126|l is valid regardless the value of K (in particular, in 
the limit for K <C 1 we have an extra small parameter, 
then Ajj(u) ~ 1 and Eq. (HI can be further simpli- 
fied). This is the power of the resonant approximation. 

We can show that the same result is achievable start- 
ing with Eq. H28f) an d neglecting the gradient terms in 
Eq. (|28|l (including, therefore, the entire y-polarization 
contribution) and also the constrained particle motion in 
that part of the phase <&t which follows from the Green's 
function, that can be found in the second term in Eq. 
(|29|l . In fact, with these prescriptions, Eq. (|112|l and Eq. 
(|113[l reduce to 



Introducing the detuning parameter C 
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Eq. l|130f) can be written as 
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Note that Ej_ has maximal magnitude when C + 
oj9 2 /2c = because otherwise the integrand in Eq. I|132|) 
starts displaying oscillatory behavior. For 9 = this 
simply means C = 0. Condition 



(133) 



is called the resonant condition. 

With the help of Eq. (fTT^I . Eq. (jH2jl can be trans- 
formed to 
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and 
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The only non-zero terms are for m = and m = —1 so 
that 



(128) 
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where the argument of Ajj is implied. Eq. (|135fl can be 
(129) integrated leading to 



Eq. (|128fl can be rewritten as 



E 



ueK 



— exp 



N m X m /2 



2 7 c 2 z J- Nw \ w/2 
( 10 



dz' < exp 



exp 



K 2 



r z ' i 

~2c~ ~ ~(3c 



' 8 7 2 fc, 



■ sin(2fc„,^') I + $ 



x 



E\ 



KojeL w e i9 » sm[CL w /2 + ujL w 9 2 /(4c)] 
c 2 ^ CL w /2 + wL w 9 2 /(4c) 



Ajjx 



(136) 



Eq. I|136(l is equivalent to Eq. ((126(1 . This can be seen 
noting that the argument in the resonant term can be 
written as 
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as it can be readily shown. The previous result can be 
cast in a more compact way, suitable for further manip- 
ulations, using the following normalized quantities 
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c 2 z jE ± 
KuieL w Ajj 
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z=— . (138) 

As remarked before, introduction of similarity techniques 
involves recognition of natural quantities which enter in 
the normalization of equations in one possible way only 
and fit with the physical nature of the problem. For 
instance, after the introduction of the normalized an- 
gle 9 we can say that 9 is to be compared, naturally, 
with the angle [c/fuL^)] 1 / 2 and not with 9 C . When 
9 ~ [c/^Lu,)] 1 / 2 , Eq. i(l"2T|) gives, naturally, the nu- 
merical factor in front of 1/N W that we were not able to 
specify before. In particular, substitution in Eq. 1127|) 
gives the following accuracy for our calculations: 
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Using Eq. (fT3"8) . Eq. fT35|) and Eq. fTTTHf) can be written 
respectively as 



E\ 



/•1/2 




/ dz! exp 


;H)- 


1-1/2 





(140) 



and 
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which are valid in the limit for z ^>> 1, 4ttN w 3> 1, and 
C + 9 2 /2 <C 4ttN w . The last relation comes from the 
fact that we are interested in the fundamental harmonic 
of our device; then Awi/oji <C 1; use of Eq. (|137(l gives 
then C + 9 2 /2 < 4wN w . Note that for C + 9 2 /2 < 1, 
that is near resonance, the formation length of the system 
becomes Lf ~ L w , as one can readily see imposing that 
the phase of the integrand in Eq. 1(141(1 be of order unity. 

As already remarked, although Eq. I|14(J|) an d Eq. 
(11411 are well known in literature (at least in their di- 
mensional form) , our discussion is far from being a mere 
repetition of textbook material. In fact we put particu- 
lar attention to the assumptions used in order to obtain 



them, which are completely neglected in textbooks. Our 
method helped to clarify these assumptions. Now, after 
this step, Eq. \ 140(1 and Eq. 1(141(1 constitute the starting 
point for further manipulations: the investigation of the 
near-field effects treated in the next Paragraph, and the 
derivation of the field in case of an electron with offset 
and deflection, treated in Paragraph IIV Dl 



C. Near- field effects 

As it has just been said, Eq. I(141|) is derived under 
several assumptions: z a ^> 1, 4wN w ^> 1, and C+9 2 /2 <C 
4ttN w . Here we will relax the condition z Q 3> 1 and we 
will treat near field effects, showing how we can control 
the accuracy of calculations using our method. 

We start considering the near field region, where z D > 
1 /2. In practice, the field at the very end of the undulator 
is not interesting because detectors are, in practice, never 
put at z — 1/2. Therefore in practical situation we will 
be always interested in situation where z Q > 1 or z Q > 2. 
In this region, suitable for our investigations, the assump- 
tion 4itN w 3> 1, and the resonant approximation are still 
retained valid, and allow substantial simplification of the 
equation for the field at the fundamental harmonic, much 
like Eq. I(141|) is a simplification of Eq. I(122|) . Using the 
same line of reasoning in the last Paragraph, but without 
expanding expressions (z — z')^ 1 around z Q , we obtain 
the following simplified expression for the field, that is 
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xexp 
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to 



2{z a 



(142) 



Here we introduced normalized units x a = 
x [uj/(L w c)} 1/2 and y Q = y a [uj/(L w c)} 1/2 . These 
definitions for normalized units, together with the ones 
given in ((138(1 , are naturally dictated by the system itself, 
through the non-normalized equations which describe 
it. In particular, x a and y a can be derived from the 
definitions of 9 and z in 1(1381) . Using normalized units 
allows one to compare each physical quantity with its 
natural measure and has the advantage of reducing the 
number of parameters that the system depends on to the 
few fitting the physical characteristics of the problem. 
Note that because of the resonant approximation, the 
phase factor in Eq. ((142(1 is much smaller than 4irN w . 

Eq. ((142(1 is valid for any value of i with an accuracy 
l/(47r7V,„). We may wish to push analytical investiga- 
tions further: in fact we have now full control over the 
expansion of (z Q — £') , meaning that we can decide 
when to truncate the series 
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both in the first factor and in the exponential factor of 

Eq. <nm 

Neglecting terms of order higher than n = m in the 
expansion of the first factor in Eq. I|142f> can be done 



with an accuracy better than (i'/i ) m+1 
that, with this accuracy, we have 



m+l 



SO 



1/2 



„=0 Z ° J -l/2 



xexp 



w2 



2(z - f) 
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Now we should study the exponential factor. Let us keep 
orders up to I = j : 
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where 9 2 = x 2 /z 2 + y 2 /z 2 . The first term to be ne- 
glect will give a contribution to the integrand equal to 
exp[i§ 2 z' j+1 /(2zi)]. Now, if (z )- j ■ 9 2 /2 < 1, this expo- 
nential contribution can be expanded too, and the domi- 
nant term in the expansion, after unity, will be just of or- 
der 9 2 /{2z ). Neglecting this term can therefore be done 
with an accuracy (z ) ■ 7 ■ 9 2 /2. If we impose that this 
accuracy be of order l/(z ) m+1 , we have automatically 
that (io)" 3 • 2 /2 < 1 and we find the useful condition 



j—m—l 



(146) 



This can be used to retain important terms: for instance, 
given a maximal observation angle of interest and an ac- 
curacy to be reached, for a certain setup, we find j such 
that (z ) ■ 2 /2 is the desired accuracy; then we find 
to such that Eq. (|146fl is satisfied. To give a numeri- 
cal example, if we choose z ~ 5 and 9 2 /2 ~ 1 and we 
want to get our result with an accuracy of about 4% we 
can put j — 2; of course we should compare this 4% ac- 
curacy with the accuracy of the resonant approximation 
l/(AirN w ), which is usually about 1% or smaller. Then, 
solving Eq. (|146fl for m we find m = 1 . Note that expan- 
sions up to m = 1 and j — 2 are easy to solve analytically. 
One obtains the following result: 



Ei =e l 



x \ [ di'exp [i (Uz' + Wz' 2 )] 

[J-l/2 

1 r 1/2 - - 1 

— / dz'z' exp [i {Uz 1 + Wz' 2 )] \ (147) 

z o J-l/2 J 



where 



U = C 



(148) 



and 



9 2 

W = — 

2z a 



(149) 



The parameter W is closely related to the near field pa- 
rameter W = L 2 w 9 2 l{2\z ) introduced in In f act ' 
once translated in dimensional units, W = 2-kW . We are 
therefore able to reproduce the results in 0], but this 
time accounting, thanks to our approach, for accuracy 
and applicability region of the approximations. 

The integrals in Eq. (|147fl can be calculated analyti- 
cally leading to 
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where 
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and 
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(152) 



It might be remarked that the erf function can be rep- 
resented in terms of Fresnel functions which were used 
in presentation of results in [44|. For the sake of com- 
parison with 0| it may be interesting to plot | B\ \ 2 as 
a function of 9. When terms in B2 are negligible this 
represents, in normalized units, the intensity in the near 
field. This is the case, for instance, if we choose z ~ 5, 
9 2 /2 ~ 1 and we want to get our result with an accuracy 
of about 20%: then we may set j = 2 and m = 0. Results 
are plotted in Fig. (|13f) and in Fig. I|14f> as intensity as 
a function of U for different values of W. 

The understanding of the region of applicability in- 
volved in the derivation of Eq. (|150f) allows us to make 
interesting remarks. For instance, note that for any given 
negative value of the detuning parameter C there is a 
value of 9 such that U = 0. This means that even when 
C = —40, for example, the value 9 2 /2 = 40 is such that 
U = and the contribution of the near field W is dom- 
inant. In this case, even choosing z Q = 40, W — 1 and 
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FIG. 13: Intensity, in arbitrary units, from the term Bi, as 
function of U for W/(2n) = 8 (upper plot) and for W/(2tv) - 
4 (lower plot). 
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FIG. 14: Intensity, in arbitrary units, from the term B\, as a 
function of U for W/(2n) = 2 (upper plot) and for W /(2tt) = 
(lower plot). 



letting m = and j = 2, the term in is the only 
contribution up to an accuracy of 5% (note that this is 
consistent with Eq. 1)146(1 . which gives a maximal angle 
9 2 /2 ~ 40. This is a very particular region of parameters, 
where the electric field depends on the near field param- 
eter W even though z D 3> 1. In this particular situation 
we can write 



E\ 



-e«. 3^ + erf 



(153) 



D. Electron motion with offset and deflection 



The same remark given in Paragraph lM Bl for electrons 
on a circular trajectory is valid here for particles in an 
undulator: the nominal trajectory is just an approxima- 
tion. Electron beams have always some finite geomet- 
rical emittance and they can be thought, in agreement 
with paraxial treatment, as a composition of perfectly 



collimatcd beams with different deflection angles with 
respect to the orbital plane of the nominal trajectory. 
As remarked in Paragraph IIII Bl this representation will 
be of great importance in the calculation of complicated 
quantities like the field autocorrelation function, which 
is of uttermost importance in the characterization of the 
statistical properties of a light source. Let us now dis- 
cuss how to calculate E from a single particle moving 
in an undulator with a given angular deflection and off- 
set with respect to the orbital plane of a nominal elec- 
tron. Once we answer this question for a single particle, 
we can add up contributions from different electrons the 
way we want, in perfect symmetry with Paragraph lHI Bl 
Let us introduce electron deflection angles r\ x and rj y and 
displacements l x and l y . The particle velocity will be 
characterized, as a function of z' as 



v ± (z') 
so that 
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X + [VyVz] V , (154) 
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one has 



if 2 
'8 7 2 fc„ 



A. 



% + % 



few 7 



sin (2k w z') H — ^cos(fc tu z') . (157) 



We will now work under the approximations: z Q 3> 1, 
4:irN w 1 and we are interested in the fundamental 
harmonic of our device. We can use a procedure analo- 
gous to the one used in Paragraph IIV Bl fand Paragraph 
IIV Cp to get the following simplified expression for the 
field in normalized units in resonant approximation: 



where 





(158) 



-K (159) 



The same normalization for 8 holds also for Tlx,y 
The accuracy of Eq. (|1 58|) is given by the ratio of 

ix/Z ~7) X ^j + 1/2 (Oy - ly/Zo - f/y^j tO 



C+l 



/2 (6 



4nN w . When C < 1 and (^ jy — l x . v /z — f] XiV ) 2 < 1 the 
accuracy is just ~ 1/(4^^). Finally, Eq. (|158|l can be 
integrated giving 



In the limit when l x/y /z <C 1, one can simplify further 
Eq. thus getting: 



where 



(el + §l) Zo , 
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and 
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Vy 



(162) 



(163) 



(164) 



(165) 



The same remarks given for Eq. H79|l hold for Eq. (|160|l . 
and we will repeat them here, since we consider them 
of great importance. Eq. I|160f) is an extremely useful 
tool, because it describes the radiation from an electron 
with offset and deflection as in an electron beam with 
finite emittance, including the correct phase factor for 
the field. Eq. Q79J1 was derived here from first principles. 
An alternative derivation based on the intuitive picture 
that the undulator radiation can be approximated as a 
sum of spherical waves emitted at the entrance of each 
pole is given in |47j . Accounting for correct phase means 
that contributions from different electrons can be sim- 
ply added up to give the field from a beam with given 
emittance at any observer position. Starting from Eq. 
(I160JI then, it is possible to calculate the field correlation 
function, and to provide a study of transverse coherence 
properties of the radiation from a given electron beam by 
means of analytical techniques. Again, in perfect anal- 

a 



ogy with what has been remarked in Paragraph IIII Bl 
numerical code can always be developed, either starting 
from Eq. I|28|) or just from the Lienard-Wiechert fields, 
which calculates the field correlation function in a generic 
case, but such a code would not help in physical under- 
standing of the situation. On the contrary, Eq. 1)1601) 
includes all relevant information about an electron in a 
realistic beam (i.e. with offset and deflection) and, be- 
ing an analytically manageable expression, constitutes a 
first step towards the more ambitious goal of characteri- 
zation of transverse coherence properties from undulator 
radiation, which we leave for future work. 



where $y is given in Eq. I|159|) . while 




(160) 
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E. The effects of electron beam emittance on the 
basic characteristics of undulator radiation 

In the next Section we will make use of Eq. I|160|) 
to calculate radiation from a complex system: we will 
demonstrate in this way the power and the practical con- 
venience of our computational method. In this Para- 
graph instead, we will study in detail how to apply Eq. 
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(|160fl to calculate angular distribution and frequency 
spectrum of radiation produced by an electron beam in 
a standard undulator when electron beam emittance is 
present. We restrict our attention to the asymptotic case 
of a large horizontal and a small vertical emittance. This 
limiting situation is, in fact, of great practical interest 
for today's third generation light sources and it will be 
discussed here in order to illustrate the effectiveness of 
similarity techniques. 

First we assume that we can use Eq. ((162J) in place 
of Eq. (|16U[I . Indicating with (5 ox , y the minimal values 
of the betatron function in the x and y directions and 
with h x , v (r]x,y) the angular distributions of the particles, 
again with respect to horizontal and vertical directions 
we have 



Intensity 




h x , y (Vx,y) = J^ P exp [ - lx 2 v ] , (166) 

V27TOV, B ' y Za x >. y > J 

where cr 2 ,/ , = £ x ,y/Pox,y Upon introduction of 



FIG. 15: Solid line : relative intensity I from a beam in the 
limit Ex > oo } Ey ► 0. The intensity is plotted as a function 
of C + Oy/2. Dashed line : single particle intensity I as a 
function of C + §1/2 for X = 



fio — L w Po , 
£ = (w/c)£ , 



(167) 



where we neglected all indexes in x,y,x' and y' for nota- 
tional simplicity, Eq. i|166|) can be rewritten as a function 
of normalized quantities as 
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We will treat the limiting case e y j f} oy <C 1. When j3 oy ~ 
1 as it usually is the case, this is simply a condition on 
the normalized omittance e 9 < 1. 

In this particular case, the particle density distribution 
in the vertical phase space behaves like a ^-distribution. 
Then, the beam intensity is simply given by 



N p c 



h = 

4iry / 27ra x i 
or, in normalized units 



drj x \E±\ 2 e-xp 



o 

Tx 
2al 



(169) 



h = 



where 



1 



2-Ke x /(3 , 



I — oo 




j ~ sin 2 C 
ar\ x 



V x Pox 

24 



h 



c rz j 



4tt 



N p c \KueL w Ajj 



h 



(170) 
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Eq. i|170|) simplifies further in the limit e x /f3 ox ^> 1. 
Also in this case, when (3 ox ~ 1 as it is usually verified, 
this is a condition on the normalized emittance i% ^> 1. 
In this limiting situation, the function sin 2 £/£ 2 in Eq. 
(|170|l behaves as the (5-distribution 8(9 x —t) x ) when com- 
pared with the exponential function, so that Eq. \ 170(1 
can be simplified as 



exp -6 2 x (3 ox /(2e x ) ^ sin 2 c 
h = L , J / df, x ^± (172) 



Now we can introduce the relative intensity 
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so that 
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and according to our approximations 
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With an obvious change of integration variable f] x 
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FIG. 16: Solid line : approximation of / using Eq. 11781 . 
Dashed line : relative intensity I from a beam in the limit 
£ x — ► oo, £ y — * 0. The intensity is plotted as a function of 
C + §1/2, as the solid line in Fig. [TSl 



FIG. 18: Relative intensity I from a beam in the limit £ x — > 
oo, £ y — > 0. The intensity is plotted as a function of 6 y for 
different values of C. The symbol C in the figure legend 
should be read as C > 0. 
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FIG. 17: Solid line : relative intensity I from a beam in the 
limit £ x — > oo, £ y —* 0. The intensity is plotted as a function 
of By for (7 = 0. Dashed line : single particle intensity J as a 
function of 0y for 9 X = and (7 = 



FIG. 19: Relative intensity I from a beam in the limit £ x — > 
oo, £ y — > 0. The intensity is plotted as a function of 9 y for 
different values of C < 0. The symbol C in the figure legend 
should be read as C. 



fj x we finally obtain 



sin 
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(176) 



It is now evident that similarity techniques we used are 
not only a theoretical tool: in fact I is a universal func- 
tion of (7+0^/2 and Eq. (|176fl can be integrated with the 



help of simple numerical techniques to get the universal 
plot in Fig. E| (solid line). This universal plot is com- 
pared, always in Fig. ^1 (dashed line), with the relative 
intensity in the case when also i x /$ox < 1. In this case 
we simply have the single particle intensity 



sin 



I = 



6/2 + eyA + 91/4, 



(177) 



C/2 + 02/4 + ^ /4 
The dashed line in Fig. ^| refers to the case 9 X = 0; 
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note the different normalization factor between I and /. 
It may be useful to provide a fit of the universal plot in 
Fig. El usm g the following function: 



Aexp 



B 



Z > Co , Co < 



£ < & < o 



(178) 

Fitting for A, -B, £i, £2 and £3 and using, for instance, 
a reasonable value for £0 = —3.00 we find A ~ 5.7, 
S ~ 0.24, a ~ 1.7, £2 ^ 2.8 and £3 ~ 1.4 • 10 3 which re- 
sult in the solid curve in Fig. El the dashed curve being 
a comparison with the real universal function plotted in 
Fig. El It should be noted that has no particu- 

lar physical sense: it simply gives an analytical approx- 
imation for the universal plot in Fig. 1151 with graphical 
accuracy shown in Fig. 1161 

After Fig. El is tabulated, one can plot the intensity 
as a function of 9 y for given values of C simply solving 
C + 8y/2 = X for 0y, X being any of the tabulated 

abscissas with constraint X — C > 0. By fixing C = 
we obtain the universal plot in Fig. El (solid line). This 
universal plot is compared, always in Fig. 1171 (dashed 
line), with the normalized intensity as a function of 9 y 
in the case when also e x /(3 ox -C 1. The dashed line in 
Fig. 1 171 refers to the case 9 X = 0; again, note the different 
normalization factor between I and /. Note that, for this 
case, Eq. (|177ll calculated for C = already constitutes 
a relatively good approximation for the plot in Fig. 1171 
Letting C vary we obtain universal plots of the inten- 
sity parameterized with respect to C; we show some of 
them in Fig. 1181 and Fig. El with respect to C > and 
C< 0. 

Finally, with the help of a Monte Carlo simulation de- 
scribed in more detail in the following Section^ we com- 
pared our results with the case of a beam with different 
values of i x and i y , always assuming /3 ox — f3 oy = 1.0 
and no energy spread is present. Note that in general, 
the spectrum (at 9 XtV — 0) or the angular intensity (at 
C = 0, 9 X = 0) is a function of two parameters: e x / fi ox 
and ty/p o y. Results are shown in Fig. I2TJI and Fig. |^ 
where relative intensities are plotted as a function of C 
for 9 X — and 9 y — and as a function of 9 y for C = 
and 9 X — 0, respectively. Fig. 1201 and Fig. 1211 illustrate 
the accuracy of the asymptotic limit Eq. (|176fl . 



V. SET OF TWO UNDULATORS WITH A 
FOCUSING TRIPLET IN BETWEEN 

In this Section we use our understanding of approxima- 
tions and region of applicability discussed before to deal 
with a difficult setup constituted by two or more undula- 
tor segments separated by strong focusing quadrupoles. 
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FIG. 20: Plots of relative intensities J as a function of C for 



and y = 0. Here /3 ox = fi oy = 1.0. 



Curve 1 : relative 

-» 00, i y — > 0, as 



intensity I from a beam in the limit i x — > 00, iy 
in Fig. 1151 Curve 2 : relative intensity from a beam with 
£ x = 40, i y = 1. Curve 3 : relative intensity from a beam 
with i x — 40, i y —> 0. Curve 4 : relative intensity from a 
beam with i x — 80, i y — > 0. Curve 5 : relative intensity from 
a beam with i x = 80, i y = 1. 
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Plots of relative intensities / as a function of 9 y for 
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and C 



0. Here (3 ox = (3 a 



1.0. Curve 1 : relative 
iy — > 0, as 



DC'. 



intensity I from a beam in the limit i x 
in Fig. 1171 Curve 2 : relative intensity from a beam with 
i x = 40, iy = 1. Curve 3 : relative intensity from a beam 



with i x — 40, £ 
beam with i x — 
a beam with i x 



y — > 0. Curve 4 : relative intensity from a 
80, i y — > 0. Curve 5 : relative intensity from 
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FIG. 22: Undulator system geometry: two undulator seg- 
ments are separated by strong focusing quadrupoles. 



In this Section, and in Fig. 1221 wc consider the particular 
case of two undulator segments separated by a focusing 
triplet. 

This scheme is known as segmented undulator scheme. 

Insertion devices at PETRA III are due to work ac- 
cording to an o?i top injection scheme and must therefore 
be compatible with injection of new particles during op- 
eration. In general, injection operations require a larger 
beam acceptance with respect to an already circulating 
beam and this poses a strict limit to a safe choice of the 
length of the insertion devices. 

Since no focusing element is present inside the undu- 
lator the (3t function increases quadratically with the 
distance from the minimal value in the middle of the un- 
dulator. In practice, then, the maximal length of the in- 
sertion device is also fixed and one has the rule of thumb 
L w ~ Pt- Here we can safely talk about the vertical di- 
rection only, because the vertical (3t function is also the 
shortest. 

In principle, given a certain acceptance £ aC c, it is pos- 
sible to increase the /^-function by increasing the un- 
dulator gap according to /3t ~ °~gap £ acc- However, the 
particle energy and the radiation wavelength are nearly 
fixed at PETRA III. Moreover restrictions related with 
the number of photons required and with operation of 
third harmonic scheme apply, which allow only values of 
K around unity. Then, \ w is fixed as well as the mag- 
netic field on axis and, as a result, the undulator gap too. 
This means that an upper limit on the @t function and 
therefore to the insertion device length L w is also fixed. 

In practical situation, when enough experience with 
the machine, safety margins might be bent. For in- 
stance a single 25 m-long in vacuum undulator oper- 
ates at SPring-8 in on top injection mode with a vertical 
Pt — 10 m. However, safety limits indicate, for a first 
design, much more stringent constraints. 

For a first conservative design, the segmented undu- 
lator scheme constitutes a reliable alternative for filling 
the first straight section at the new arc of PETRA III, 
which can accommodate a 20 m long insertion device. 
Breaking a long insertion device in several parts allows 
one to " keep the average vertical /3^-function rather small 
which in turn is essential to allow for a small geometri- 
cal vertical aperture and magnetic gap" (cited from > 



Section 4.1.8). "The maximal vertical /3-function along 
this undulator will be about 10 m. Therefore a vertical 
inner size of 9.5 mm is necessary for the corresponding 
chambers" (cited from [ll|, Section 3.2.1). 

With this motivation in mind we start our investiga- 
tion. We are interested in computing the radiation in- 
tensity seen by an observer at angles 9 X and 8 y . We will 
be interested in both angular and spectral intensity dis- 
tribution. Although it is possible to take advantage of 
the analysis performed in Paragraph llV C| and treat near 
field effect, in this paper we will restrict our attention to 
the far field limit. Then we will compare our results with 
the case of an uniform undulator. 

Of course radiation intensity can be obtained, numeri- 
cally, from the Lienard-Wiechert expressions for the elec- 
tromagnetic field without approximations of any kind. 
The previous statement though, sounds like declaring 
that solution to all electromagnetic problems can be 
found solving Maxwell's equation: it is undoubtedly true, 
but quite generic. The most general expressions for the 
Lienard-Wiechert field can be used, and have been used 
(again, see |2jj and |28j for instance) as a basis for numer- 
ical codes, but understanding of correct approximations 
and their region of applicability can simplify many tasks 
a lot, including practical and non-trivial ones like the one 
we are going to study here. Instead of relying on compu- 
tational power we propose a simple technique for dealing 
with the segmented undulator scheme which minimizes 
numerical difficulties and computation time and can be 
implemented straightforwardly also by non-expert pro- 
grammers. In order to reach this result we will use our 
understanding of approximations in undulator radiation 
theory, thus simplifying equations as much as possible 
before using, in the final step, some numerical integra- 
tion technique. It is worthwhile to stress the fact that 
we will have full control of our approximations at every 
step, meaning that we will be able to get the radiation 
intensity and the accuracy with which this is computed 
as well. 

We will consider the undulators tuned at the funda- 
mental frequency. Therefore wc start reminding the res- 
onant approximation seen in Section llVI In particular, 
using our method we were able to demonstrate that under 
the resonant approximation the field contribution from 
non-resonant structures can be neglected with respect to 
the undulator field with an accuracy of 1/ (4irN w ) . Here 
we will restrict our attention to a region of parameters 
where such an approximation is valid. Having indicated, 
as in Fig. [221 with L s the straight section length, wc 
will therefore be interested in the fundamental harmonic 
of the system within the assumptions: z Q 1 + L s /L w , 
AttN w > 1, and C + 6 2 /2 « 4ttN w . 

Also, for the same reason and with the same accuracy 
we can neglect the constrained motion in the Green's 
function phase and consider the particle moving on a 
straight line when calculating this phase term. 

In Paragraph llV Dl we treated the case of a single elec- 
tron with offset and deflection. If one is interested in 
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the far field limit case, one can use Eq. (|162|l to obtain 
the contributions from each segment but special atten- 
tion must be taken in summing up the fields with the 
correct relative phase factor. 

With the help of Eq. I|162() and the coordinate sys- 
tem depicted in Fig. |22we start accounting for the field 
contributions from the first and the second undulator in 
the case there is no straight section between the two seg- 
ments that is L s , defined as L s /L v 
we have: 



is zero. In this case 



B, - 1 r i(i'a+i>o) r -»Ci/2 sm Ci/ 2 
1 2 6 Ci/2 



(179) 



technology (see j45|); in the case of PETRA III, they 
will be usually designed and tuned in order to provide 
matching condition for particles at nominal energy mov- 
ing at zero angle with respect to the undulator axis, and 
with zero detuning. This means that the phase shifter 
will contribute for a relative phase <fi s h such that 



2c 7 2 



= 2nn 



(186) 



where n is an integer number. Taking out the unessential 
2ixn contribution (or setting n — 0) and accounting for 
tfish one gets 



while 



where 



E„ - 1 r »(f 3 +$o) r tf2/2 Sm C2/2 

2 "2 6 C2/2 
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and 



u>iL a 



2c7 2 <x>i 



Vxlf + g 



%i) 
(187) 
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Vyl 



and 



C2 = Y + \ - + ^ (0y - % 2 ) 
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(182) 
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Here 4> s is defined as usual by Eq. (|159[) : offsets I 



l c(x,y) 

refer to the middle point of the system (i.e. the point 
with zero value of z in Fig. I22JI and deflections T)(x, y )i 
and f)( Xt y)2 refer to the first and the second undulator 
respectively, and the relation between the two is deter- 
mined by the focusing elements only: therefore it makes 
sense to talk about deflections for the first and the second 
undulator without specifying the point. 

The next step is to account for the phase shift in the 
case L s ^ 0. On the one hand, shifting the first segment 
of a quantity —L s /2 gives a phase contribution: 



loL s 



2cj 2 



(184) 

On the other hand, shifting the second segment of a quan- 
tity L s /2 leads to an extra phase term: 



A$ 2 



- U!l 



2c7 2 <x>i 



2 \Py - % 2 J 
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Finally, we make use of the relation C = Alu/(2^c) and 
we express A$i and A$2 in normalized units thus writ- 
ing 



A$i 



and 



C 



t: 



Vyl 



(189) 



A$ 2 = %f. 



where 



Vx2 



7 [Oy-Vifl 



L' 



l + K 2 /2 



(190) 



(191) 



We can account for a given energy deviation A7/7 ex- 
panding C around the nominal energy. Trivial calcula- 
tions show that, for u> = u)\, u>i being now calculated at 
8 = and at nominal energy we obtain a shift Ae in the 
detuning parameter given by 



A$ 9 



2c 7 2 



Vx2) + ~ (Oy - rj y2 y 



(185) 

To complicate the situation further, we should account 
for the fact that, in order to practically control the rel- 
ative phase between the two contributions phase shifters 
will be installed. These devices are also useful in FEL 



A7 



(192) 



Including Ae in our equations we obtain the final results: 



I p i(* s +*o) p -iCi/2 p iA* 1 Sm Cl/ 2 

2 6 Ci/2 



(193) 
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and 



E'2 



_ 1 *(<£>s+*o) p <2/ 



2 tA<s> 2 sinC2/2 
C 2 /2 



(194) 



TABLE I: Set of parameters for the segmented undulator 
scheme 



with 



Quantity 



Unit 



Value 



C2 = 



A$ 2 = 



H 4- An- 1 /-> \2 1 / „ \ 2 

-^-^ + 4(^-^1) +4(^-^1) .(!95) 
2-^ + 4(^-^3) +4(^-^2) ,(196) 



C 



-L>- 



4 V 



Tjxl 



4 V 



Vyl 



(197) 



C + A E 



L'- 



Vx2j +-T[0y-Vv2 



(198) 

where C is now the detuning parameter for a particle 
with nominal energy and other quantities are as defined 
before. 

As we have said, focusing elements simply change the 
electron trajectory, but they do not contribute directly 
to the total field at the observer position. We will treat 
them like a thin lens with focal length f x and f y in 
the horizontal and vertical planes respectively. We can, 
therefore, introduce normalized focal lengths f x — f x /L w 
and f y = f y /L w . 

From now on, when possible, we will neglect indexes 
x or y for notational simplicity. If we have offset and 
deflection l\ and 771 at the entrance of the system (where 
z = —L2/2 — 1/2, see Fig. we can easily calculate 



N w 

L s 

7 
A 

for* 

u 

fry 

Ey 
fy 

C>E 

K 



m 
m 



20.0 
690 
6.5 
11742.0 
2.5 ■ 10" 10 
0.218 
24.005 
0.044 
1.0 ■ 10~ 9 
5.0 
1.291 
10.272 
0.260 

1.0 ■ 10~ n 
5.8 

1.1 ■ 10~ 3 
1.66 



Here N p is the number of particles in the 
for notational simplicity, we are taking 
as implicit. Upon introduction of 



(200) 

beam and again, 
indexes x and y 




(199) 



It is worth to note here that if we take the limit for 
f — ► 00 and L s — > 0, the sum E\ + E2 gives back Eq. 
(|160|l . as it must. 

Since we have many particles here we can consider a 
collective description of the beam, for instance in terms 
of the Twiss parameters ax, Pt and 7t for the horizontal 
and vertical planes at the entrance of the system, where 
7t = (1 + aj.)/f3r- Knowing the Twiss parameters and 
the emittance e at the entrance of the first undulator it 
is easy to write down the number of particles with offset 
between li and li + Sli and with deflection angles between 
r]i and 771 + Sr]i , always at the entrance of the system; in 
fact, one can simply write it as q(U , m }8U 8m , where g 
is the following density distribution 46]: 



It = L w j T , 
£ = (uj/c)s 



(201) 



we can write Eq. (|200|l as a function of dimensionless 
quantities as 



g(h,fji) = ex P 



7t?i + 2arhm + Prfjl 

2i 



(202) 



The local energy spread of the beam will be assumed to 
be a gaussian function with standard deviation &e'- 



l{h,Vi) = 7, — ex P 

27T£ 



7t?i + 2a T /i?7i + (3 T rjl 
2e 



,-A%/(2a%) 



y/27 



(203) 
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TABLE II: Set of normalized parameters for the segmented 
undulator scheme 



&Tx 


0.218 


PTx 


1.20 


JTx 


0.873 


ix 


25.1 


fx 


0.250 


&Ty 


1.29 


Pt v 


0.514 




5.19 


e v 


0.251 


fy 


0.290 


Is 


0.325 


O 'E 


9.54 


Li 


0.137 



where <je = 4kN w <je and 




brackets indicating the ensamble average over the particle 
distribution. 

A set of realistic normalized parameters for the seg- 
mented undulator is given in Table |U and Table [H] for 
dimensional and dimensionless quantities, respectively. 
From these tables we can already foresee a large effect 
of the normalized energy spread parameter &e which is, 
in our case, comparable with the normalized horizontal 
omittance. 

We can now average the normalized intensity over the 
particle distribution function and the energy spread of 
the beam at the entrance of the system thus finding the 
final expression for the total intensity 



I b = \E 1 +E 2 \ 2 = 



df] x 



dA E F(A E ) / dl xl 

) J — oo 

dlyi I dfjyi g x {i x i,fi x i)g y (lyi,fi y i) 



e ^Ci/ V A *! sin + e ifo/2 e iA* 2 sin C2/2 



Ci/2 



C2/2 



As in Eq. (|171|l , the normalized intensity is related to 
the non-normalized intensity lb by 



h = 



4-7T / C ZoJ 



N p c \KojeL w A 



jj 



We now reached an expression, Eq. I|205|) , which can be 
integrated with the help of numerical techniques. Results 
of numerical integration are the intensity as a function of 
the detuning parameter C for fixed observation angles or 
the intensity as a function of angles for a fixed detuning 
parameter. 

We developed a short code which uses a Monte Carlo 
method to compute Eq. (|205|1 . A quasi-random gener- 
ator is used to fill the horizontal and vertical transverse 
phase space with a given number of particles. The ac- 
tual simulation makes use of Sobol sequences to reduce 
the formation of cluster of particles, a well known phe- 
nomenon described, for instance, in 48] and references 
therein. Then, the intensity for each particle is calcu- 
lated and all contributions are summed up to give an 
estimation of Eq. (|205[1 . The critical parameter in this 
kind of codes is the number of particles used. In general 
the higher the number of particles is, the smaller the fluc- 
tuation of the results from run to run and the higher the 
accuracy of the integration. On the other hand numerical 
algorithm designers warn that quasi-random sequences 
start to display unwanted correlations after several tens 
of millions calls to the generating routine. We observed 
that a choice of several million particles is more than suf- 
ficient to ensure that results converge fast enough so that 
(204) no relevant fluctuation is present in the outcome (with 
graphical accuracy of the results which will be described 
in the next few figures). 

First, we compared results from the Monte Carlo ap- 
proach in the limiting case L s — > 0, f x . y — > 00, e x — > 00, 
i y — > and oe — ► 0. Results were in perfect agreement 
with the ones shown in Fig. El Fig. El Fig- El an d 
Fig. E| as it should be. Then, as has already been said 
in Paragraph IIVDI we compared this asymptotic situa- 
tion with a more practical case of a uniform undulator 
with e x ~ 40 and i y ~ 1. Results have been already 
described, and shown in Fig. 1201 and Fig. EHl 

Finally we used all the features of our program in or- 
der to deal with the segmented undulator case. We set 
up parameters as in Table [I] which correspond to a de- 
sign set for the PETRA III case. The solid line in Fig. 
1231 represents the intensity as a function of the detuning 
parameter C at X = and 9 y = 0. The dashed line in 
Fig. 1231 is the intensity for the uniform undulator case, 
corresponding to the limit L s -*0m and / — > 00, while 
L w = 20 m. In this case we set the minimal betatron 
function values to f3 ox — 20 m in the horizontal direction 
(205) and j3 oy = 10 m in the vertical one. Beam energy, wave- 
length, undulator deflection parameter, relative energy 
spread and electron beam emittances are left unvaried 
as in Tabled Assuming the minimal (3t function in the 
center of the undulator we can calculate the Twiss pa- 
rameters at the system entrance using for both horizontal 
and vertical plane (see, for instance, 0|) : 



(206) 



Moreover t) 2 is expressed in terms of l\ and 7)1 in Eq. 
(E!J. 



Pt{z) = Pa 



0o 



(207) 
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FIG. 23: Solid line : relative intensity from a beam with 
parameters listed in Table [n] The intensity is plotted as a 
function of C for 8 X — and 9 y = 0. It is normalized to 
the maximal intensity from a beam in a uniform undulator. 
Dashed line : relative intensity from a beam in a uniform 
undulator (L s — > m and / — * oo). We set the minimal 
betatron function values f3 ox = 1 in the horizontal direction 
and j3 y = 0.5 m in the vertical one. i x , y and &e are left 
unvaried as in Table UTI 
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FIG. 24: Solid line : relative intensity from a beam with 
parameters listed in Table [n] The intensity is plotted as a 
function of 9 y for 9 X = and (7 = 0. It is normalized to 
the maximal intensity from a beam in a uniform undulator. 
Dashed line : comparison with the relative intensity from 
a beam in a uniform undulator (L s — * m and / — * oo). 
We set the minimal betatron function values f3 ox = 1 in the 
horizontal direction and f3 oy — 0.5 m in the vertical one. i x , y 
and &e are left unvaried as in Table ITT1 



1t{z) = 7c = ^ , (208) 

a T (z) = -7 (2 - z w ) , (209) 

where z w indicates the z-position at the beam waist, 
while subscripts " o" , as usual, denote minimal values of 
the Twiss parameters. Taking z w = m at the center of 
the system and z = —10 m we obtain the initial Twiss pa- 
rameters to enter in our program: ax x = 0.5, (3tx = 25.0 
m and jtx = 0.05 m -1 , ax y = 1.0, (3t v = 20.0 m and 
7T y = 0.1 m -1 . 

It is important to remark that we normalized the solid 
line in Fig. that is the result for the segmented un- 
dulator, to the maximal intensity found for the uniform 
undulator while the dashed line represents the relative 
intensity for the uniform undulator. This allows direct 
comparison between uniform undulator and segmented 
undulator schemes. 

To conclude, in Fig. [2] (solid line) we plot again the 
intensity as a function of 9 y at C — and 9 X = 0. Simi- 
larly as for Fig. [231 the dashed line in Fig 1231 represents 
a comparison with the uniform undulator case L s = m 
and / — > co, where the other parameters are left unvaried 
as in Table HJ Normalization of the curves are as in Fig. 
12*51 The final results for the two schemes do not appear 
to be dramatically different. 

As a comment for both Fig. [231 and Fig. [21 it is 
interesting to note the smoothing action of finite beam 
omittance combined with finite energy spread. 

VI. CONCLUSIONS 

In this paper we have shown that even well-known 
and long studied subjects, like Synchrotron Radiation, 
continue to retain interest from a theoretical viewpoint. 
In particular we have proposed a technique to compute 
harmonic contributions of the electric field from a mov- 
ing charge particularly suitable for analytical investi- 
gations. We showed that our general expression, Eq. 
(J2BJ, can be obtained from already known expression 
Eq. J2J) with the use of a paraxial approximation, but 
also much more straightforwardly from Maxwell's equa- 
tion in paraxial form, with the help of a parabolic Green's 
function method. The main advantage of our approach 
with respect to standard techniques is that we separated 
from very beginning approximations independent from 
the system under study (paraxial approximation) and as- 
sumptions that, from time to time, arise in the study of 
a particular magnetic system. This gave us the possibil- 
ity to specify region of applicability and accuracy of well 
known results, which we were, of course, able to recover 
using our method, like the field from a charge moving 
on a circle, in a short bending magnet or in an undu- 
lator and also edge radiation effects from a particular 
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setup. Again, our particular approach is logically much 
less involved with respect to the others, since it relies 
on direct integration of the paraxial Maxwell's equations 
for a given harmonic. As a result, our expression gives 
a clearer and completely fresh insight in old results. Be- 
sides being able to discuss their region of applicability 
and their accuracy, we pointed out that phase corrections 
with respect to the spherical wavefront for an electron 
on a circle arise naturally from analytical treatment. We 
have seen that short magnet radiation is a complete sep- 
arate phenomenon with respect to edge radiation. We 
stressed the importance of the knowledge of the entire 
electron trajectory in order to obtain a correct calcula- 
tion of the field: from this viewpoint we have seen how 
edge radiation from two dipoles separated by a straight 
section arises, in our method, from the contribution due 
to the straight section alone, and that the field term from 
the magnets is completely negligible. We discussed, un- 
der resonance approximation, accuracy and applicability 
region of the far field undulator radiation and we dealt 
with near field effects as well: in particular we were able 
to reproduce already known results studying their region 
of applicability in detail, and we found also other regions 
of interest, where the near field parameter W is of fun- 
damental importance. 

Besides this, our method proved to be a reliable basis 
for calculation of more complicated quantities like field 
autocorrelation functions, either from bending magnets 
or undulators, which are of fundamental importance in 
the understanding of spatial coherence properties of Syn- 
chrotron Radiation. In fact, by direct application of our 
general formula Eq. (|28|l we were able to provide ana- 
lytically manageable expressions for bending magnet or 
undulator radiation in the case of an electron with a 
transverse offset and deflection. These expressions can 
be integrated over realistic beam particle distributions 
and can be used in very practical calculations. Their 
strength stems from the fact that they correctly account 
for the field phase and that they can be simply summed 
up to give results for apparently complicated situations. 



Complete treatment of spatial coherence would proba- 
bly double the size of this paper. For this reason, this 
first work focuses on the basis of Synchrotron Radiation 
theory only. 

A few applications which exploit the power of our ap- 
proach have been selected for the final part of the paper. 
At the end of Section ITVl we addressed basic character- 
istics of undulator radiation in the presence of electron 
beam emittance. We discussed the asymptotic situation 
for a large horizontal and a small vertical emittance. 

Finally, in Section we considered a novel undulator 
configuration (setup of two undulators separated by a fo- 
cusing element) which is planned for installation at PE- 
TRA IIL Computation of radiation characteristics from 
new setups by means of numerical techniques alone re- 
quires, almost always, modifications of existing simula- 
tion codes which can be better done by the code authors 
themselves. Up to date, the particular case we chose to 
study has not been included in existing simulation codes. 
For this reason we selected this particular example: our 
goal was to demonstrate that, in general, the application 
of our method allows solution of apparently difficult prob- 
lems relying on simple computer algorithm which can be 
developed by non-expert programmers, with obvious ad- 
vantages in saving time and better physical understand- 
ing. 

We expect that our technique will be used in the future 
as a basis for developments of simulation codes and for 
physical understanding of complex situation as well. 
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